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NUMERICAL VALUES OF TWO-CENTRE INTEGRALS 
FOR 3D ELECTRONS 


By MACIEJ SUFFCZYŃSKI 
Institute of Theoretical Physics, University of Warsaw, Warsaw 
(Received May 15, 1961) 


Nuclear attraction integrals and overlap integrals for 3d electrons in chromium, manganese, 
iron, cobalt, nickel, copper and zinc have been computed. The atomic wave functions have been 
approximated by analytical expression with two exponentials. 


Calculations based on the tight-binding approximation may be of some use in the 
study of 3d electrons in transition metals. The computation of the necessary integrals iıı the 
two-centre approximation has been suggested by Fletcher and Wohlfarth (1951). The calcu- 
lation of nuclear attraction integrals in crystals is different from the case of the free molecules. 
To simulate the crystal potential one calculates the nuclear attraction integrals with the 
atomic sphere excised (Fletcher 1952). The radius of the sphere is always taken equal to 
half the distance between the neighbours, for simplicity. 

We report here a numerical study of the two-centre integrals, both nuclear attraction 
and overlap, for 3d electrons in chromium, manganese, iron, cobalt, nickel, copper and 
zine (see Suffezyhski 1959). The numerical computations were done on Ferranti Mercury 
Computer. 

The atomic potential seen by a 3d electron was approximated by U(r)= Z,(7)/r, with 
Z,(r) = 1+ (Z— 1) e-?”. Here Z is the atomic number of the atom, and y for all the con- 
sidered elements was taken equal 3. This was a simplification dictated by a lack of better 
knowledge of the actual potential for the 3d electron. 

'The radial part of the atomic 3d wave function has been approximated by analytical 


expression with two exponentials only. 
P,() = (Bye + Bye?) (1) 


The coefficient B, by the smaller exponent 5, was chosen so as to fit the tail of the 
numerical wave function. The choice of b, affects the tail usually slightly, (see Suffczyński 
1955). A b, was chosen rather arbitrarily. Then the machine calculated B, so as to get an 
exact normalization of the radial wave function according to eq. (3) of the author's reference 


(1956 c). 
(045) 
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The use of the two exponentials only in eq. (1) is the most obvious drawback of the 
present caleulation. Different choices have been tried for the exponents, resulting n 
ical functions which fit in different ranges of r the numerical wave functions given in 
the literature, The use which can be made of the numerical results reported here requires 
first of all to decide which analytical wave function to choose for a problem at hand. At 
the time when this work was being done the report of Watson (1959) was not yet 
accessible. 

In the present calculation of the nuclear attraction integrals (E) the exact formulae for 
the two-centre energy integrals in the crystal have been used, as published by the author 
(1956 a, c). In the notation of ref. (1956c) all the terms v= 0, 1, 2, 3, 4, 5, 6, 7 have been 
computed and summed here. Concerning the formulae of ref. (1956c) it should be noted 
here that the (ddr) energy integral will have the proper + sign if the coordinate systems 
for the two centres will both be made right-handed. Further a misprint in the formula 
(21) of ref. (1956c) should be corrected at the present occasion: the term 27015 should 
read 2700/a5 (see Suffczyński 1959). 

The formulae for the overlap integrals for 3d electrons are published by Kotani et al. 
(1955) and by Jaffe (1953) and Roberts and Jaffe (1957). The overlap integrals (S) have to 
reduce to unity for zero distance between the two centres as we use normalized wave func- 
tions. 

The lattice constants quoted in different sources are slightly different. For that reason 
the two-centre integrals have been calculated for different distances distinguished by sub- 
scripts. The order of the neighbour is denoted by a roman numeral. The half-distance be- 
tween a neighbour and the central atom is denoted here by R and expressed in Bohr units, 
The atomic constants have been taken from the tables of Cohen and DuMond (1955). 


l Bohr unit of length — 0.529172 x 107? cm 
l atomic unit of energy — 27.20976 electron-volts. 


The energy integrals (E) are expressed in electron-volts using this atomic unit of energy. 
In the tabulation of the results the necessary letters describing the integrals are written 
once for each metal and are omitted in subsequent numerical cases. The values of the two- 


behind the comma is the power of ten by which the number before the comma should be 
multiplied. In the tabulation only the absolute values of the integrals are given, in other 
words, we list in each case 


—E(3d3de) ^ E(3d3dm)  — E(3d3dà) 
S(3d3de) —S(3d3dm) S(3d3d0) 


Chromium, Z=24 


The selfconsistent calculation of the atomic functions of 0৮ without exchange has 
been published by Mooney (1939). The 3d radial function can be approximated. with the 
two exponents ړا‎ 2, b3—5 (see Suffczyński 1955. An error in B's given there for chro- 
mium is corrected below). It was tried also to approximate the 3d function of chromium by 


واد نادد 
r YY‏ 
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an expression similar to the one approximating 3d function calculated with exchange for 
Mn** by Hartree (1955). This led to the exponents 5,—2.3, b,—5. 

| The lattice constants in the body-centered cubic chromium is quoted as a,=2.879 x 
x10 ® cm. The corresponding half-distances of the three consecutive neighbours are in 
Bohr units, 


a,V3/4 = R(L) = 2.355838 
a,j2 = RüL) = 2.7202875 
a/V2 = RIL) = 3.847067 


The lattice constant of the b. c. c. chromium is also quoted as aq=2.8845 x 10-8 cm. 
The corresponding half-distances of the neighbours are then, in Bohr units, 


2/314 = R(I,) = 2.3603387 


6212 


= REL) 


= 2.7254843 


qall/2 = R(III,) = 3.854417 


The two-centre integrals have been computed for both lattice constants and they are 


listed here to show the variation with interatomic distance. 


B, — 3.448 B, — 49.7046 
(3d3do) (3d3dm) (34340) 
1 E 7.150, —1 4.853, — 1 8.1072 
S 6.863, —2 4.840, —2 8.897, —3 
Er E 3.694, 17 1.817, —1 2516,22 
| S 4.042, —2 2.173,—2 3.211, 23 
| 
m E 2.222, —2 6.697, —3 5.956, —4 
S 3,055, —3 1.204, —3 1.093, — 4 
F E 7.682, — 1 4.796, — 1 7.992, = 2 
5 6.825, —2 4.795, —2 8.788, —3 
II, E 3.582,—1 1.791, —1 2.474, —2 
S 4.008, —2 21172 3.164, —3 
8 E 2.178,—2 6.547,—3 5.809, —4 
5 3.886,--3 1.179, —3 1.068, —4 
b, c 2 b a 5 
B, = 3.558 B, = 46.39986 
I, E 8.007, —1 5.103, —1 8.554, —2 
Û 7.249,--2 51292 9.442, —3 
Il, E 3.806, —1 1.917,—1 2.660, —2 
5 4.283,—2 2.306, —2 3.411,—3 


UC t Ut Hy ob 


SUP 


০ t Ut Ut Hy la ty اما‎ 


& bg BSF 


০০ by 


Um Bp Ut nN 


5 


= 16.76491 


5 
6 


1.50236 


2.3592,—2 
4.205, —3 
8.027, —1 
7.209, —2 
3.763, — 1 
4.247, —2 


2.306, —2 
4.132, —3 


TOTO 9 
1.048, —1 
5.279, —1 
6.333, —2 
3.456, —2 
6.369, —3 
1.070, 0 
1.043, — 1 
5.222, — 1 
6.280, —2 
3.389, —2 
6.259, —3 


3.601; —1 
3.029 —2 
1.295, — 1 
1.374, —2 
3.755, —3 
6.469, —4 
3.558, — 1 
3.002, —2 
7 
1.358, —2 


3.664, —3 
6.327, —4 


7.100, —3 
1.280, —3 


5.043, — 1 
5.081, —2 


1.890, — 1 
2.279,—2 


6.941, —3 
1.254, —3 


7.144, 1 
1.583, —2 
2.746, — 1 
3.450, — 2 
1.054, —2 
1.944, —3 
1.063, — 1 
7.514, —2 
2.707,—1 
3.409, — 2 
1.031, —2 
1.904, —3 


1.714, 1 
1.601 —2 
5.051, --2 
5.721, —3 
9.169, —4 
1.599, — 4 
1.689, — 1 
1.582, — 2 
4.962, —2 
2.635, —3 


8.924, —4 
1.559, —4 


6.320, —4: 
1.162, —4 


8.433, —2 
9.326, —3 
2.616,—2 
3.361, —3 


6.164, —4 
1.136, —4 


1.228, —1 
1.411, —2 
3.877, —2 
5.136, —3 
9.441, —4 
1.767, —4 
1.210, —1 
1.394, —2 
3.813, —2 | 
5.060, —3 ı 


9.210, —4 : 
1.727, —4 | 


2.303,—2 | 
2.348, —3 | 
5.602, —3 | 
6.846, —4 | 
6.778, —5 | 
1.204, —5: 
2.264, —21 
2.313, —3! 
5.578, —3 

6.725, —4 
6.582, —5 | 
1.172, —5: 
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Manganese, 2-25 


The atomic functions for Mn** have been calculated with exchange by Hartree (1955). 
The 3d radial function was approximated with two exponents b, = 2.3, b,— 5. 

Manganese has a complicated crystal structure. In the caleulation of the two-centre 
integrals the distance between the nearest neighbours was taken d — 2.24 x 10-8 cm. The 
corresponding half-distance is, in Bohr units, d/2— R(I)= 2.1165141. 


B, — 4.64 B, — 61.50236 
(3934০) (39347) (34399) 
T E 6.829, — 1 3.762, —1 5.756, —2 
S 4.699, —2 l 2.999, —2 5.141, —3 
Iron, 2-26 


For iron the wave functions calculated by J. H. Wood (1955) have been used. Wood 
was taking the exchange effects into account by using an average exchange potential pro- 
posed by Slater (1951). The 341 (spin up) wave function computed by Wood was approxi- 
mated by several analytical expressions with one exponent near 2 and the second near 5. 

In the b.c.c. iron the lattice constants a,= 2.8663 kX and d,= 2.4823 kX were 
measured by Basinski, Hume-Rothery and Sutton (1955). Taking 1 kX = 1.002063 x 10-8 cm 
we get the corresponding half-distances of the first order and second order neighbours, 


in Bohr units, (Du Mond and Cohen 1953) 


dj? =R(L) = 2.345471 
dı l2 = R(IL) = 2.7082937 


The lattice constant for the b. c. c. iron is usually quoted as a, = 2.86645 x 10-8 cm. 
The corresponding half-distances of the three consecutive neighbours, are, in Bohr units, 


(Cohen and DuMond 1955) 
a,/3/4 = R(l)) = 2.345569 
)اه‎ = R(II,) = 2.708429 
a, V2 = 70711) = 3.830297 


(343৫০) (3d3dn) (343৫9) 

E 4.119, —1 2.214,—1 3.328, —2 

S 31832 1.954, =2 3.289, —3 

II, E 1.688,—1 7.456, —2 9.410, —3 
5 1.654, —2 8.051, —3 1.102, —3 


SUO bj Oby اما‎ 


= 5.1 
, = 82.55929 


joe) 


= 5.1 
= 77.80095 


ماما ly‏ اما BE Dh‏ بح ما بحم ېخ 


$ Ut 


t2 
I 


5.1 
79.44393 


R 


hh Dy & bij 


Su 
ټم‎ 


=5 

B, = 68.22085 
E 

5 


4.118, 1 
3.133, —2 
1.687, —1 
1.654, —2 
1.949, —3 
1.251,—3 


2.503, —2 
1.215, —2 
3.308, — 1 
2.502, —2 
1.249, —1 
1.215,—2 
4.394, —3 
1.230, —4 


3.709, — 1 
2.877, —2 
1.416, —1 
1.408, —2 
3.708, —1 
2.877, —2 
1.416, — 1 
1.408, —2 


3.210, — 1 
2.466, —2 
1.174, —1 
1.155, —2 
3.713, —3 
6.146, —4 


3.202, —1 
2.604, —2 


2.214, —1 
1.953,—2 
7.453, —2 
8.049, —3 
2.090, —3 
3.548, —4 


1.423, —2 
5.443, —3 
1.638, — 1 
1.423, —2 
5.107, —2 
5.441, —3 
1.136, 3 
1.913, —4 


1.860,—1 
1.648, 2 
5.842, —2 
6.334, —3 
1.859, — 1 
1.648, —2 
5.840, —2 
6.332, —3 


1.543, —1 
1.350, —2 
4.654, —2 
4.990, —3 
9.321, —4 
1.575,—4 


1.543, —1 
1.375, —2 


3.327, —2 
3.288, —3 
9.406, —3 
1.101, —3 
1.733, —4 
3.026, —5 ` 


2.225, —3 
6.948, —4 
2.299, —2 
2.225, —3 
6.034, —3 
6.945, —4 
8.868, —5 
1.534, —5 


2.630, —2 
2.588, —3 
6.944, —3 
8.10154 
2.629, —2 
250019 
6.941,—3 
8.103, —4 


2.107, —2 
2.044, —3 
5.351, —3 
6.174, —4 
1.079, —5 
1.226, —5 


2.058, —2 
2.018, —3 
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II, E 1369,—1 4.517,—2 5.067,—3 
S 1179, 2 4.916,—3 5.900, —4 
WERA b, 2'51 
B, = 4.09 B, = 15.52839 
I, S 2.464,—2 130198 T9183 
115 5 Ilo, 2 4.654, —3 5.588, —4 


Cobalt, 7=27 


The 3d wave function for cobalt was taken the same as for iron. The lattice constants 
for hexagonal cobalt are a = 2.5074 x 10-8 cm and c = 4.0699 x 1075 cm. The resulting 
half-distances of the first order neighbours are, in Bohr units 


4 c(3/8)* 77077222579 

4 (a?/3 + c2|4)^ = R(1,) = 2.359667 
2/2 = 18075) = 2.369173 
2ه‎ = R(Ij) = 3.845536 


The half-distances for the second and third order neighbours are in Bohr units, 
alV2 = R(II) = 3.3505165 
a/3|2 = R(III) = 4.1035279 


In the face-centered cubic cobalt with the lattice constant dy = 3.5442 x 10-8 cm the 
half-distances of the first and second order neighbours are 


ł a,//2 = R(I) = 2.367971 
aJ2. = R(II,) = 3.348817 


b, =2.1 bę 25.1 

B, — 2.815 B, — 81.97126 
(3d3do) (3d3dn) (3d3d6) 

7 E 4.081,--1 23255] 3.244, —2 

S 3.087, —2 11011 3.198, —3 

fe E 4.035, —1 2.144,—1 3.191, —2 

S 3.064, —2 1.890, —2 3.154,—3 

a E 3.045, —1 2.084, —1 3.088, —2 

S 80157 1.848, —2 3.066, —3 

7 1 7.251,—3 1.995, —3 1.645, —4 

S 1.203, —3 3.390, —4 2.876, —5 

II E 3.028, —2 1.001, —2 9.755,—4 

S 4512753 15429,—3 1.457, —4 


b, 54 
B, — 11.80095 


Ej ما بح ما‎ bg ما‎ by ما بح ما بح ما‎ ই 


: 
eee 


C t by. 


51- 
82.55929 = ر 


CH 


৩০ 


بح ما بخ ما بح UC b ০2‏ 


a bj & bû 


3.956, — 1 
3.023, —2 
3.042, —2 
4.144, —3 


3.667, —1 
2.829, —2 
3.622, —1 
2.805, —2 
3.020, —1 
2.157,—2 
4.860, —3 
8.108, —4 
2.239, —2 
3.078, —3 
3.545, — 1 
2.763, —2 
2.250, —2 
3.002, —3 


3.272, —1 
2.460, —2 
3.232, — 1 
2.439, —2 
3.153, —1 
2.397, —2 
4.208, —3 
6.927, —4 
1.951, —2 
2.637, —3 
1.842, —3 
3.309, —4 
3.163, — 1 
2.402, —2 
1.961, —2 
2.648, —3 


2.092, —1 
1.853, —2 
1.007, —2 
1.436, —3 


1.822, —1 
1.609, —2 
1.795, —1 
1.590, —2 
1.742, —1 
1.552, — 2 
1.254, —3 
2.135,—4 
6.931, —3 
9.022, —4 
1.749, —1 
1.557,—2 
6.972, —3 
9.973, —4 


1.606, — 1 
1.389, —2 
1.582, — 1 
1.373, —2 
1.535, — 1 
1.340, —2 
1.081, —3 
1.822, —4 
6.007, —3 
8.484, —4 
4.350, —4 
1.942, —5 
1.541, —1 
1.344, —2 
6.042, —3 
8.528, —2 


3.101, —2 
3.077, —3 
9.814, —4 
1.465,—4 


2.559, —2 
2.512, —3 
2.515, —2 
2.475, —3 
2.429, —2 
2.402, —3 
9,756, —5 
1.704, —5 
6.366, —4 
9.493, —5 
2.440, —2 
2.411, —3 
6.407, —4 
9.549, —5 


2.230, 22 
2.160, —3 
2.200, —2 
2.128, —3 
2.124, —3 
2.065, —3 
8.393, —5 
1.453, —5 
5.497, —4 
8.108, —5 
3.124, —5 
5.827,—6 
2.134, —2 
2.073,—3 
5.532, —4 

j 8.156,—5 
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3,=23 b, =5 
D B, — 68.22085 
E 3.248, —1 1.510,—1 1.200,—2 
S 2.556, —2 1.840,—2 1.05778 
E 3.206,—1 1.486, —1 1.964, —2 
S 2:53] 22 1.323, —2 1.926, —3 
E 30199074 1.439,—1 1.893, —2 
5 2.483, —2 1.289, —2 1.866, —3 
E 3.198, —3 7.758,—4 5.711,—5 
5 5.355, —4 1.323, —4 9.959, —6 
E 1627 =2 4.728, —3 4.104, —4 
5 2.952523 6.781, —4 6.108, —5 
E 1.332, —3 20714 2.024, —5 
5 2.428, —4 5.479, —5 3.798, —6 
E 3.133, —1 /1445,—1 1.902, —2 
S 2.489, — 2 1.293,—2 1.873,—3 
II, E 1.636, —2 4.757,—3 4.131, —4 
5 2.262, —3 6.818, —4 6.146, —5 
E bE b, = 5.1 
B, = 4.09 B, = 15.52839 
Ti 5 DALE oo? 1.268, —2 1.852, —3 
80 5 2.395, —2 125212 1.823, 23 
p S 2.349, —2 1220,22 1.766, —3 
iP 5 5.074, —4 1234 9.440, —6 
II, S ae 6.423, —4 5.788, —5 
T S 2.355,—2 1:2247522 EÊ 
Ii, S 2.143, —3 6.459, —4 5.823, —5 


Nickel, Z— 28 
For the 3d function in nickel the approximation assumed by Fletcher (1952) has been 


used. 
The lattice constant of the f. c. c. nickel is quoted as a,= 3.52394 x 107? cm. The 


corresponding half-distances for the first and second order neighbours are, in Bohr units, 


2 به‎ ۷2 = R(L) = 2.35443 
2112 = RI) = 3.3296735 


If the lattice constant were taken as a5 — 3.608 x 10-8 cm the half-distances would be 


2 azj/2 = R) —241 
a,2. = R(II,) = 3.408254675 
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R(L) is the half-distance used in the original calculation of Fletcher (1952). 


D b, =5 
B, — 1.979 B, — 85.86026 
(3d3do) (3437) (3৫3৫9) 
is E 3.787, —1 2112-7 3.302, — 2 
5 2.651, —2 1.765, —2 3.151, —3 
Il, E 3.424, —2 1.208, —2 1.250, —3 
5 4.400,--3 1.654, —3 1.813, —4 
i E 3.339, —1 1.805,—1 2.745,—2 
II, E 2.780, —2 9.511, —3 9.573, —4 


The appoximation of Fletcher consists in taking only the terms with v= 0, 1, 2, 3 in 6 
integration over all space and only the term with v= 0 in the integration over the atomic 
sphere, in the energy integrals (see Suffezynski 1956c). This approximation gives in the 
present case the values 


1; E 3.3316, — 1 1.8006, — 1 2.7392, —2 


which are very close to the values obtained from the full expressions with y= 0, 1, ...7. 
Thus the approximation of Fletcher, that is omitting the integrals with unequal exponents 
over the atomic sphere, has been found, in the course of the present computation, to be very 
satisfactory indeed. 


Copper, Z= 29 


The wave functions for atomie Cut have been calculated with exchange by Hartree 
(1936). It has been approximated by Fletcher (1952) and the same, essentially, approxima- 
tion was used here. 

The lattice constant of copper is a= 3.608 x 1073 cm. The corresponding distances 
of the first and second order neighbours are, in Bohr units, 


łaj/2 = R(I) = 2.4] 
৫12 = R(II)= 3.40825 


b, =2 15125 
324.979 D, = 85.86026 
(3d3do) (3d3dz) (3d3d6) 
T E 3.383, —1 1.823, 1 2.766, —2 
S 2.443, —2 1.562, —2 2.695, —3 
TL E 2.800, —2 9.562, —3 9.614, —4 
S 3.705, —3 1.343, —3 1.427, —4 


Zinc, Z = 30 


The wave function for 3d electrons in zinc has been taken the same as in copper. 
The lattice constants for hexagonal zinc are a = 2.659 x 10-8 cm and c= 4.935 x 10-8 cm. 
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The resulting half-distances for the first order neighbours are, in Bohr units, 
aj2 = R(L)= 2.5124156 
4(a?/3 + c?/4)? = R(I,) = 2.7459173 


b = 55 
B, = 1.979 B, = 85.86026 
(3৫3৫০) (3d3dn) (3d3d6) 
I, E 2202, 2]. 1.374, —1 1.981, —2 
S 2.085, —2 2022 20772 
I, E 1.544, —1 6.986, —2 9.067, —3 
S 1.404, —2 7.192, —3 1.030, —3 


In the present calculations the numerical values of Fletcher (1952) and of the author (1956b) 
have been satisfactorily confirmed. This gives confidence in the accuracy of the former and 
in the programme of the present numerical work. 

My thanks are due to. Mr. S. Michaelson for introducing me to Ferranti Autocode 
programming and to Miss R. Fenn for her effective help in checking the programmes. 
The work was done in 1959 in the Department of Mathematics, Imperial College, London. 
The access to the Ferranti Mercury Computer of the London University Computing Unit 
is gratefully acknowledged. 
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THE SPINOR SPACE AS AN EUCLIDEAN COMPLEX SPACE 


By JERZY LUKIERSKI 
Institute of Theoretical Physics, University of Wrocław, Wrocław 


(Received May 23, 1961) 


The equivalence between Rzewuski's spinor space and 4-dimensional Euclidean complex 
space is proved. The free field equation in spinor space is solved by the procedure 01 analytical 
continuation. The field equation with sources is investigated, using the complex distribution Ö(z). 


Introduction 


In a series of papers Rzewuski (1958a, b, 1959, 1960a, b) investigated the spinor space 
as geometrical basis for the theory of elementary particles. The full group of transformations 
in spinor space is a direct product of two 6-parametric unimodular groups C x C’ and 
a phase factor group a. The simplest faithful representation 22) (a, P = 1, 2) of the 12-para- 
metric group C x C' was used as the complex coordinate of spinor space. 

Recently, Mozrzymas and Rzewuski (1961) proposed an interesting method of solving 
a second order differential equation in spinor space. This equation, in the paper by Moz- ' 
rzymas and Rzewuski, is treated as an analytical continuation of the Klein-Gordon equation, 
and its solution in spinor space can be obtained from any solution of the Klein-Gordon 
equation by the procedure of analytical continuation. In the present paper, this method of 
solution is shown to be intimately related to the isomorphism between the spinor space and 
Euclidean four-dimensional complex space. 

The equivalence between spinor and vector spaces was investigated by Rzewuski (1958b) 
and, independently, by Roman (1960). Rzewuski introduced an eight-dimensional real 
vector space with two real metric forms, isomorphie to the original spinor space. In this 
paper, we consider complex vector spaces. It is shown that the interpretation of spinor 
space as an Euclidean complex vector space admits of explaining certain results, hitherto 
calculated explicitly, as simple consequences of this property of Rzewuski's space-time theory, 


l. Isomorphism Between the Spinor Space and the Euclidean Complex Space 
Let us introduce a complex four-dimensional space with Euclidean metric: 
z, = 2? = inv. (1.1) 


The linear transformations 
(up = 1, 2,3, 4) (1.2) 
(957) 


958 


are limited by virtue of (1.1) as follows: 


aa, = 0, (dm — complex). (1.3) 


pv up 


For Rzewuski’s spinor space 230 we have 


লাক =u? - inv (০১19, y, 6 = 1,2) (1.4) 
where 
لاو‎ NE. 152 = لوو‎ কির رر په وزو‎ 
z. = fgg وهن تي په‎ ই ۳ Min জা ال‎ 
If 
Û yj, CIT 
Zap Sz 9 B (1.5) 
we have, by (1.4), 
i 5%] = Sêr = | (1.6) 


The space with the invariant forms (1.1) and (1.4) are connected as follows: 


„aß (1.7a)‏ ښوه 


Zu £ 01193 
OT 

28১8 = CAB (1.7b) 
where the 0;,,g are the elements of the well-known Pauli matrices, 044,5 = i054. The 
justification of (1.7) results from (1.1), (1.4). After substitution of (1.7a) in (1.1) or (1.7b) 
in (1.4), we have 


bo 
ro 


u? =z? (1.8) 
The isomorphism between the spinor and complex Euclidean spaces may be stated for 


infinitesimal transformations. From (1.2), we have 


62, = daS, (1.9) 
and, from (1.5), 
Óz,,8 = 092 Zy,g + Sq; 05 (1.10) 
The relations (1.3) and (1.6) imply 


0৫, = — 0a,, (1.11) 
and 
061 =— 958) 6ق‎ 1 = — 552 (1.12) 
By (1.7b), we have 
65৫)8 = 0//2)992% (1.13) 


Using the properties of ede we obtain directly from (1.9-10) and (1.13) that 
day = 4 (6 a ge 8021 75798 „052 js (1.14) 


We see that the infinitesimal parameters of transformations in spinor and complex 
Euclidean spaces are connected linearly. 
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On assuming 
0a,,— (da,,)*, (1.15) 


- we obtain an Euclidean space with complex vector components and real transformation 
coefficients. From (1.14-15), we obtain 


05 ——0S*t 6S'=—6S't (1.162) 
where ÒS, ö5” are infinitesimal 2x2 matrices (see (1.10)) which implies for finite 
Sos Pug Sidus] (1.16b) 


The Minkowski case (complex vectors, six-parametric group of Lorentz rotations) is 


determined by the choice 


da; = (da,)* day = — (da,)* (1.17) 
Using (1.14) and (1.17), it is readily shown that 
OS 09° (1.18a) 
which is fulfilled only if 
9 . (1.18b) 


The cases for which 6S=0 or 05'=0 were investigated by Rzewuski (1958b). Stat- 


ing the results of Rzewuski’s paper in our notation, we have 


ö5“ —0: day = ep 9৫74 (1.19a) 


ij 

085 = 0: day = — eg, da (1.19b) 

Thus it is shown that the Euclidean complex space includes spinors in C (and scalars 

in C’) as well as spinors in C’ (and scalars in €). By a special choice of the infinitesimal 


transformations in spinor space we can obtain the result that the spinor coordinates 2,5 
describe the components of the complex Euclidean or Minkowski vectors. 


2. Solution of the Free Field Equation in Spinor Space 


Let us investigate only second order differential equations in spinor space. We have the 
following two second order differential operators, invariant with respect to the group C x €': 


9 2: TA = (2.1) 
D = =r en RB 
There are two points of view on the theory of interaction in spinor space (Rzewuski 


1960a): 
1) The free field equation (together with the initial conditions) determines the free 


fıeld uniquely. 
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In this conventional case, the differential equation includes D as well as D*. An example 
of this type can be provided by the following equation, justified by the symmetry proper- 


ties: 
{F(D, D*) + F(D*, D)} f(z „aza; g) = 0 (2.2) 


where F is a polynomial in D, D*. 
2) The free field equation admits the degrees of freedom describing interaction. 


E.g. the equation investigated by Rzewuski and Mozrzymas (Mozrzymas 1961, Rze- 


wuski 1959, Rzewuski 1960a, Mozrzymas and Rzewuski 1961) 
(70687425828) fms as0 (2.3) 


admits of arbitrariness in the dependence of f (39, 226) on Ze 
In this paper will be continued the latter treatment of the free field, which represents 
a new approach to the problem of interaction. We assume the following free field equation: 
SMS SCNT Els : 
(D—x* M (26৩ )} (5858, 2858) (2.4) 
where M? is an arbitrary scalar constructed. from "BP 
By (1.7a) we have, instead of (2.4), 


9? bU V. 
LR des cn Fun. SJ "U (2.5) 
The substitution 
Z,= M(ziz, (2.6) 
implies 
E 
lazê = al Fia Lp) = 0 (2.7) 
CLu 
if ZZ is complex. 
In the case of 
Mîz dan (2.8) 
as realized in eq. (2.3) (4 = An, we have 
az, (2.9) 
and it is necessary to introduce a new real variable & 
* = 
gs a. =) (2.10) 
Zn 


as investigated by Rzewuski (1959). Because of 


TREE . em) 
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where 
9Z, *9 Oa xm 
1 14022) টি (2.12) 
we have using (2.6) 
9 3 A 2 
= eel zay — al. 
টে | PAREZE (2.13) 
and eq. (2.5) can be rewritten as follows: 
2 PESO" | 
(sz 2 Z) 17292 (2.14) 


where Z, describes only seven degrees of freedom (see (2.9)) and the remaining degree is 
described by the angular variable a. 

The case described by the equality (2.8) (where A=) was investigated recently. 
From our considerations it is easy to deduce the method of solving the equation presented 
in the paper by Mozrzymas and Rzewuski (1961). 

Let us introduce the following vector components y, 


Yu-% 62752) (2.15) 
where 
7 oa )2.16( 
It is readily shown that 
ae (2.17) 


The condition (2.16) is invariant under the transformations in spinor space satisfying 
the relation (1.18b). Assuming (1.18b), we see from (2.17) and Sect. 1 that the y, are usual 


Minkowski vectors. 
We now write the Klein-Gordon equation in the variables y,: 


92 
[Z - el roa =0 (2.18) 


The general solution of eq. (2.4) with the value of M determined by (2.8) can be ob- 
tained from (2.18) by means of the following two steps: 


9 0 Py, © 

< -- real 2.19. 
1) 37 ৯ 3y. (m 7 « — real . (2.19a) 
A By om Se Z, complex, Z2 real (2.19b) 


The first step can be interpreted as the introduction of a certain interaction (Rzewuski 
1960a, b). The general solution of the Klein-Gordon equation after the interchange l) was 
already obtained in Rzewuski's paper (1959). The second step consists in the analytical 
continuation of the coordinates. Using (2.6), (2.8) and (2.10), we obtain the solution as 


a function of 22১, 
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The first step is connected only with the form (2.8) of M2. If M? is not bilinear in the 
length of 5: the general solution can be obtained by the procedure of analytical continuation 


Yu > Zu Zu complex, Z2 complex. (2.20) 


3. Curved Orthogonal Basis in Spinor Space 


Let us consider the Euclidean complex vector space z, (with real or complex 22). The 


transformations (1.2) can be generalized as fellows: 
> A yy (23) (3.1) 
wherein 


Aufl) Ara) — Öko (3.2) 


* . . . . . 
As the dependence of A,, on zı is not, in general, determined by a function of in- 
. * . . 
variant from z?', the transformation properties of u,; 


u, = Au) 2, (3.3) 


differ from the transformation properties of z,. The two most interesting cases are these: 
1) u„is a scalar with respect to C’ and a vector component with respect to C, and inversely 
2) u, describes a scalar with respect to C and a vector component with respect to C’. Assum- 
ing that i 


Ii NA . Û 
iy B: >: ad + SQ 5 
pm (zę) = 3 ow 3 Op; Ova; y €;ı নল (3.4) 
VE; azi 
where €,, is a unit vector in space, described by € or 
Le u żę 
zr ox E 3; rê ` e 
E pe (zo) dkg of; on Ova; y Oof; ô Cr; Vel, zê (3.5) 
A €4; 24 


where e,, is a unit vector in space-time, we obtain these two cases, respectively. Indeed, 
using the properties of g,— matrices, the transformation properties of A,, and Au are 
readily determined. We can also prove that the relation (3.2) is satisfied for (3.4) and (3.5). 

The vectors z, are defined by means of the formula (1.7a). Let us insert in (3.3) the 
vectors Z, defined by (2.6). It is easily proved that 


U, Fo» "ee Cn Z, (3.6) 
and 
| U, لال‎ Z, (3.7) 
satisfy the relations 
WRZE 08) 
U; = US: U مس‎ A 
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Comparing the bilinear formula of Rzewuski (19588, 1959) for the position four-vector 
in space-time. 


AD E A yi 
Xu; = $05 75,25, ûn Ex; (3.9) 


with eq. (3.6), and on using (1.7), we obtain the result that the vectors X „and U, are identi- 
cal. As the formula (3.6) defmes four real variables, and Z, describes seven degrees of 
freedom (assuming (2.8)), we can in addition to (3.6) define three more variables په‎ inde- 
pendent of X,, The orthogonal character of the transformation (3.6) implies 


"pue Pee a 9 9 
928 ES 9X2. I D (ke. و07‎ IK A (3.10) 


and, if the dependence of the field function on Z, can be expressed by means of X, or if 


RE =0 (3.11) 

ALA i : 
we have 

DF — 0. (3.12) 


The equality (3.10) is a consequence of the result that the variables (3.9) present a special 
case of the variables (3.3). Finally, we obtain the interpretation of X, as coordinates in 
a certain orthogonal curvilinear coordinate system defined in a real four-dimensional sub- 
space of the full complex four-dimensional Euclidean space z, (with eight real degrees of 
freedom). 


4. Remarks on the Equation with Sources in Spinor Space 
Eq. (2.4) may be written with the delta-function on the right hand side: 


Ê SW M a) a (Zu, Zu) = 0(24) 0(2) (z3) Ö(za) (4.1) 


92 
The delta-function with complex argument z can be defined as 


ö(z) = > Jk e?z dp (4.2) 
jê 


where the path of integration is chosen as follows: 
LE DZIW RER <00 (4.3) 
and the direction of the path is determined by the condition Re dp > 0. Thus, we have 


Rega) ess 
dp = Re] z* de (4.4) 


if 
da >0. 
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From (4.2), we obtain 
+0 
Re(z*) 23 1 [^ D 
= বব — | 65৫ 4.5 
ö(z) দল 2a. 61:21 da (4.5) 


and, fınally, 


Re(z*) ick a9 
Rez) ^ Ol) 


It is easy to prove that, for an arbitrary path of integration K, we have 
1 x SFE) if zọ belongs to K 
: s PAULUS lo if z, lies outside K 


K 
l 1 if z belongs to K 
ö(z — zy) dz = 
K 


ó(z) = 


(4.7) 
and 


0 if zą lies outside K 


Formula (4.6) can also be proved by help of the definition of the delta function depend 
ing on a complex argument, as introduced by Nakanishi (1958) (see Appendix). 

The analytical continuation of the delta function can therefore be effected, as, for the 
real argument x we have 


x 
O(x) = bd x ó(x?) (4.8) 
" 

Formula (4.8) is a special case of (4.6), if z = z*. We obtain the interchange ó(x)-ó(z) 
if, in the relation with real v, y’ 


Re y’ 
Ô KN EE af Ó FEU : 
(y) Rey] 7009) (4.9) 
we put 
pzy = yy (4.10a) 
Passing to complex arguments, as follows (see Fig. la, 1b): 
ya, MER! (4.10b) 
we obtain from (4.8) the equality (4.6). 
im y ım y' 


@ 


szyna — Rey 


Fig. la Fig. lb 
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By (4.6), eq. (4.1) can be written in real variables, with delta functions of real argu- 
ments. Let us introduce, for each complex component z, the variables as follows: 


gre’? (4.11) 
The differentiation operator 22 m be expressed thus: 
zZ 
= = e-ip uh Û 2 4.1 
Po or  r0g9 ae 


eit dm 2i 9° 1 ০১ 4.13‏ — ت 
f ar? r 9197  r$8g? তি‏ 924,7 

The complex delta function (4.6) can be expressed as follows: 
ö(z) = 0(g) e(r) (4.14) 


where 


lif 0<9 <, «< 9< 2x 
0 (g) = x : )4.15( 
uf 2 <? z” 


The equation (4.1), written in the variables r,, p, where! 


w W: == re"? (4.16) 


eu 


is of the form 


x Ir, Ar JG Piz dg? 


১ 
১০2৮6 ZUM Wr 9»; x?) ZĘ [I 0 (a) Ò (ra) (4.17) 


u=1 


where 


I = DO (4.17a) 
WE 


Eq. (4.1) or (4.17) defines Green's function for the equation with sources in the full 
Euclidean complex space z,, defined by X, (physical part of spinor space), and by a, and œ 


(isotopic part of spinor space). 


1 Einstein summation convention is not used. 
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We obtain the physical interpretation of the solution of the inhomogeneous equation 
in spinor space (for example (4.1)) in accordance with Rzewuski's prodecure (1960) on 
averaging over the four degrees of freedom a,, « independent of X,. The second possibility 
of a physical interpretation can be related to the average over the four degrees of freedom 
@ (see (4.16)). In this method, the solution of (4.17) can be interpreted physically after the 


following” procedure : 
1 2 
1077727222৮ Ga fa... 1 Ar, 943203) WÎ (ZE) (4.18) 
৪ 9 


where 7;, ir, define four space-time degrees of freedom. 

The interpretation in (4.18) implies that the 4-vector of position in space-time is in- 
variant only with respect to the transformations of the 6-parameter subgroup C x € of 
C xC (see (1.18b)). This property is easy to prove, taking into account the relations (2.15) 


and (4.16). 
The author wishes to thank Professor Rzewuski and Mr Mozrzymas for their helpful 


and valuable discussions. 


APPENDIX 


Nakanishi (1958) introduced by means of the Cauchy formula the following definition of 


Ók(z) = 2 (2 _ ij (A.1) 


21 


where z+} defines the path of integration K,, and z —the path of integration K_ (see Fig. 2a, 
2b, 2c). i 


© 


Fig. 2a Fig. 2b Fig. 2c 


Since 


fe) = x m 


(A.2) 


where C, is an infinitesimal circle about the point z —0, satisfying the relation C. =K,- K_ 
the properties (4.7) of the complex delta function are satisfied. 
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Let us investigate only the straight lines K. If z=0 belongs to K, the argument for an 
arbitrary point of the path is equal to gy or gy +r, where po determines K. 
We have in this case: 


T 1 = 
Da, ced AŻ a 
zo |z| ie te 
i = gi” 1 TÊ a T pa 3 
Zt re (A.3) 
and 
Rez ال‎ 1 1 
9 2) = 190 =, 
5 |Re z| 32 Qi n Xie j= >| (A.4) 
Since 
নতি a 
dni ||z|+ie — |z|— ie] * (A.5) 
for the straight line gg = const, gg + zt = const 
we have 
Reiz), 
ó(z = eP Â 
(2) = [Re(z)| ° (2) (A.6) 
Using 


902) = |z] ٠ (|z|?) 
we obtain from (A.6) that 
Re (z*) z*5(lz|2 
Formula (A.7) is valid for an arbitrary path of integration with a continuous derivative 
in the point z=0. Indeed, from (A.7) we see that only for z=0 must the form of the curve K 
be taken into account. If in (A.6) we introduce the following interchange 


Po" > (0) (A.8) 


where the formula z= |z| ۶۱م‎ describes points of the path of integration K in the infinites- 
imal neighbourhood of the point z=0, the relations (A.6-7) are also satisfied. 


ó(z) = (A.7) 
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ZUR MÓGLICHKEIT DER ELEKTRONENHAFTSTELLENANALYSE 
MIT HILFE DER EXOELEKTRONENEMISSION 


Von B. SUJAK 
Institut fiir Experimentalphysik der Universitiit, Wrocław 
( Eingegangen am 29 Mai 1961) 


Die Móglichkeit der Elektronenhaftstellenanalyse auf Grund der Erscheinungen von 
Exoelektronenemission wird besprochen. 

Es werden besonders die Móglichkeiten einer Analyse mit Hilfe der Messungen von: 
1) Postelektronenemission (Abklingkurven), 2) Photostimulierten Coelectronenemission und 
3) Thermostimulierten Coelektronenemission (Glow-Kurven) zusammenfassend diskuttiert. 


Die Untersuchungen der Exoelektronenemission sind vom Anfang an mit dem Ziele 
gefiihrt worden, diese Erscheinuag in der Zukunft zur Elektronenhaftstellenanalyse der 
dünnen Oberflüchenschichten des Festkórpers auszunutzen. 

Im Grunde genommen liegen bis jetzt drei Móglichkeiten vor, die Haftstellenanalyse ' 
einer Oberflichenshicht mit Hilfe der Exoelektronen durchzuführen, und zwar durch: 

1) Eine Analyse der Abklingkurven (Postelektronenemission). 

2) Eine Analyse des Emissionsstromes der Exoelektronen, der wihrend einer Er- 
würmung der Probe mit konstanter Geschwindigkeit der Temperaturerhóhung gemessen 
wird (thermostimulierte Coelektronenemission, ófters Glow-Kurven Analyse genannt). 

3) Eine Analyse der Abhüngigkeit der Exoelektronenemissionsintensitit von der 
Lichtwellenlinge des beleuchtenden Lichtes auf gleiche auffallende Lichtenergie (oder 
Quantenzahl) reduziert (photostimulierte Coelektronenemission). 

Diese Grundmóglichkeiten werden auf den nachstehenden Seiten umrissen. Es werden 
auch von anderen Autoren abgeleitete Beziehungen angeführt. Auf die Ableitung selbst 


wird aber verzichtet und nur auf Originalarbeiten verwiesen. 


Postelektronenemission 


Nassenstein (1955) hat die Gleichungen der Theorie der Phosphoreszenz auf die 
Exoelektronenemission übertragen, um mit diesen, die mit der Zeit in Dunkelheit ab- 
klingende Emission von Exoelektronen beschreiben zu versuchen. 

Bei der Voraussetzung, dass es sich um einen monomolekularen Prozess handelte, 


schrieb er: 
(969) 


970 


a) für den Fall der einzelnen diskreten Haftstellenarten und bei der Annahme, dass 
jedes, nahe an der Oberfläche befreite Elektron den Kristall verlassen kann 


٤ ৪ 
I= Ange RT ° exp [ —ste zr | (1) 


oder auch 


12255 (9) 
b) für den Fall der gleichmissigen Verteilung der Haftstellen, 

I=N.(kT[t) (1— چ‎ *) (3) 
wobei für Zeiten £>1077 sec (welche eingentlich in Betracht kommen) dieser Ausdruck 
praktisch in 

I-NET-17 (4) 
übergeht. 


c) für den Fall einer Haftstellenverteilung, die einem Exponentialgesetz N, = 4e^ ** 
folgt; wird der Emissionsstrom folgend beschrieben: 


I = f (s, kT) R لسا‎ 0 (5) 
dass heisst, die Emission verlüuft nach dem Gesetz: 
ټم‎ my (6) 
wobei 
a=akT+1 (7) 


Wenn die besetzten Elektronenhaftstellen in einem kleinen Energieintervall (eı, 4) liegen 
und in diesem Energieintervall N, konstant ist, verläuft die Elektronenemissionsintesität 


wie folgt: 
exp E ste | — exp | sie | 


Die Bedeutung der Symbole bei den Abgeleiteten Abhängigkeiten ist folgend: 
e — energetische Tiefe der besetzen Haftstellen in eV (thermische Aktivierungs- 


1N, EL (& 


` 


energie), 
T — absolute Temperatur, 
t — Zeit in sec, ; X 


k — Boltzmann-Konstante =8,62 * 10-5 eV/Grad 


S” exp E A — Wahrscheinlichkeit, dass ein Elektron bei der Temperatur T aus 
einer Haftstelle der Tiefe e thermisch befreit wird, 

s — etwa 108 sec}, 

Q — eine Grösse die angibt wieviele der befreiten Elektronen den Kristall als Exo- 
elektronenstrom verlassen. können, bei (2) bis (8) gleich 1 gesetzt, 

ny — Anzahl der Elektronen, die in den für die Emission wirksamen Haffstellen der 
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Tiefe e zur Zeit t=0 gelagert waren, 
B — Konstante. 
Zur Illustration des vorgesagten wird eine theoretisch errechnete Abklingkurve der 
Exoelektronenemission fiir eine konstante Haftstellenverteilung (N, = konstant) in dem 


J(0) 


NE 


10° 


107 10? 10? 


t [sec] 


Fig. 1. Theoretisch errechnete Abklingkurve der Postelektronenemission für eine konstante Haftstellenvertei- 
lung (Ng = konstant) im Energieintervall von (sı, €) &= 0,55 eV bis e= 0,75 eV aber N,=0 für €< &, 
€ > €, und für T= 300?K (nach Nassenstein 1955). 


Energieintervall e, = 0,55 eV bis &;= 0,75 eV und für T= 300°K nach Nassenstein (1955) 


in Fig. 1 wiedergegeben. 


Thermostimulierte Coelektronenemission (Glow-Kurven) 


Die Glow-Kurven-Methode beruht auf der Messung der Intensität der Exoelektro- 
nenemission während die emittierende Probe mit konstanter Geschwindigkeit erwärmt wird. 
Bei dieser Methode sind — je nach der Anzahl der Arten von mit Elektronen besetzten 
Haftstellen — bestimmte Maxima auf der Exoelektronen-Intensitütskurve zu finden. Wenn 
diese Maxima scharf voneinander getrennt zu finden sind, ist die Identifizierung und die 
Analyse der entsprechenden Haftstellen erleichtert, oft aber nur in diesem Falle möglich. 

Die Haftstellenarten werden durch deren energetische Tiefe & charakterisiert (Fig. 2). 
Durch diese wird im Falle der Glow-Kurven-Methode die thermische Energie gemeint, 
die zum Befreien des festhaftenden Elektrons aus der Haftstelle benótigt wird. Diese Energie 
wird oft als thermische Aktivierungsenergie bezeichnet. 

Zur Bestimmung der energetischen Tiefe e, der als Exoelektronen-Emissionszentren 
dienenden Haftstellen, wird eine anniherende Beziehung zwischen dieser und der Tempe- 
ratur T,„, die der Temperaturlage des Maximum entspricht, verwendet. Diese Beziehung 
wurde zuerst vom Randall und Wilkins (1945a, b) zur Bestimmung der energetischen 
Tiefe der Lumineszenzzentren abgeleitet. Sie wurde aber zuerst vom Bohun (1954) auf die 
Exoelektronen Glow-Kurven übertragen (Fig. 2). Diese Beziehung lautet: 


ET, (9)‏ 25 حم 
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oder auch, 
Tm 


e = zg [eV] (10) 


Bohun und Dolejsi zeigten auch (1959), dass die an Alkalihaloiden erzielten Messer- 
gebnisse auf eine Monomolekularen Prozess hinweisen und dass die Austrittsarbeit nur 
einen sehr kleinen Einfluss auf die Temperatur T,, hat. Die Austrittsarbeit beeinflusst ledig- 
lich sehr stark nur die Intensität der Exoelektronenemission. Es scheint, dass die Be- 


N/t 


— ae ---— Jope 9 
€, ০০ ৫7০০ 
P Tm1 
Û 500 
£5 

2 7772 

500 

a) T/500 6) 


Fig. 2. a) Energetisches Schema von 2 Haftstellenarten bei der Oberflüche eines Festkórpers. P, und P, — zwei 

verschiedene (energetisch genommen) Haftstellen, Ej — gleich Null gesetztes Energieniveau des Elektrons 

im Unendlichen (wenn in Ruhezustand), Eg — Verlauf der Energiebarriere für das Elektron an der Oberfläche : 

des Festkórpers, e£, und €, — Aktivierungsenergie der Haftstellen P, und P,. b) Glow-Kurven (schematisch) 
die dem Fall der Haftstellen P, und P, entspricht. 


ziehung (9) oder (10) nur auf das Volumen einer Wirksamen Oberflüchenschicht bezogen 
werden kann, wie die allerletzten Untersuchungsergebnise z. B. von Seidl deuten (1960). . 

Bei der Diskussion des Problemes führte Seidl (1960) zwei Interpretationsmechanisme : 
ein. Einer von denen ist mehr als eine Oberfliichenemission zu betrachten. Bei diesem | 
Mechanismus ist das Entweichen der Elektronen aus den Haftstellen lediglich durch die : 
Emission bedigt. Bei diesem Mechanismus wird e einmal aus dem Anlaufteil des Maxi- 
mum des Emissionsstromes 


N çê 
— = Ryl e M (11) 


und das anderemal aus der Richardsonschen Gerade, auf die schon bereits Bohun hinge- 
wiesen hat, 


T BRR | (12) 


ady 
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bestimmt. Dabei bedeuten: 


XS Impulsgeschwindigkeit (es wird angenommen, dass ein jedes emittierte Elektron 


registriert wird, was bis jetzt nicht nachgewiesen wurde, besonders bei den 
Messungen mit Hilfe des Luftspitzenziihlers), 


ny — Anfängliche Anzahl der lokalisierten Elektronen, 


v — Schwingungsfrequenz des Elektrons in der Potentialmulde des Haftstellen- 
zentrums, 

Q — Geometrische Wahrscheinlichkeit der Emission, 

k - Boltzmannsche Konstante, 


T — Temperatur der Probe, 

b — Geschwindigkeit der Temperaturerhóhung der Probe, 

T, — Temperatur, die der Lage des Maximum auf der Glow-Kurve entspricht. 

Die Übereinstimmung der Werte ৪, die aus (11) und (12) ermittelt werden, soll zugleich. 
ein Kriterium für die Verwendbarkeit der gegebenen Interpretation sein, das heisst, dass 
das Entweichen der Elektronen aus den Haftstellen lediglich durch die Emission bedingt 


wird. In diesem Falle kann aber auch e aus der folgenden Beziehung annäherend ermittelt 
werden: 


LB an) kim (13) 
Tp E 
wobei 7, die Temperatur ist, die der halben Höhe des Maximums auf der Seite der niedri- 
geren Temperatur entspricht. ; 
Da die Lage des Maximums T, auf der Glow-Kurve von der Geschwindigkeit der 
Temperaturerhöhung abhängt, wird sie von den Sachbearbeiter auf etwa 1 Grad/sec bis 
2 Grad/sec bei allen Messreihen festgelegt und bei den Ergebnissen stets angegeben. 
Der maximale Emissionstrom (bei der Temperatur 7) ist nicht nur von der anfänglichen 
Anzahl der lokalisierten Elektronen ny, sondern auch von der Geschwindigkeit der Tempe- 


raturerhóhung der Probe b auf folgende Weise abhingig: 


-N 70 Eb 
En = A 14 
( t Ihe e kT}, A9. 


Der zweite Fall, der vom Seidl diskuttiert wurde, ist mehr durch einen inneren Mecha- 
nismus gesteuert, wie z. B. eine Überführung des Elektrons aus der Haftstelle in das Lei- 
tungsband. Dabei kann die Austrittsarbeit fast keine Rolle spielen. In diesem Falle gelingt 


er wiederum zu der Formel (9): 
&,— C + ET, (15) 


wobei nur dann C—l ist, wenn 


1> 


2 ৩০. 
Væ DIA (1+ e kTm (16) 


b Eo 6০০ 
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Die neu eingeführten Symbole haben folgende Bedeutung: 


e, — Aktivationsenergie des inneren Prozesses, 


রি — Austrittsarbeit der Elektronen, 
Q, — Geometrische Wahrscheinlichkeit des Überführens des Elektrons in das Leit- 
fühigkeitsband, 
ي9‎ Geometrische Wahrscheinlichkeit der Emission des Elektrons nach aussen. 
Seidl bemerkt jedoch, dass die Ungleichheit (16) wahrscheinlich ziemlich hiufig er- 
füllt ist. In diesem jedoch Falle kann e, verhältnismässig genau aus den bei verschiedenen 
b vorgenommen und im System In b, 1/T,, oder im System In bj T? verarbeiteten Messreihen 


laut der Gleichung 
EI Te 
10 ne পারা 
ron 8 (1+) (7) 
b 


& = k1,, 0 


ermittelt werden. e, dagegen kann aus der Verarbeitung des Anlaufteiles der Kurve (18) 
mit Hilfe der Richardsonschen Geraden bestimmt werden 


5০০ 


RT (18) 


= POOG E 


Aus dieser Diskussion ist zu entnehmen, dass die aus dem Anlaufteil des Maximum 
bestimmte Energie (Gleichung (11) oder (18)) in der Regel der Austrittsarbeit der Elek- 


9 
£ 
///////////// 
পি Es >C 
a) 6) 
Fig. 3. Energetische Schemen der Emissionsmechanisme wobei: L — Boden des Leitfihigkeitsbandes, Eg— 


Verlauf der Energiebarriere für das Elektron, Ej, — gleich Null gesetztes Energieniveau des Elektrons in 
Ruhezustand ist. 
a) Emission wird durch die Austrittsarbeit bedingt. b) Emission wird durch die Aktivierungsenergie des inneren 
Prozesses bedingt (nach Seidl 1960). 


tronen entspricht. Die Aktivationsenergie des inneren Prozesses entspricht dem Wert von 
£p, der mit Hilfe der Gleichungen (9), (10), (12), (13), (15) oder (17) ermittelt wird. 

Für die Na- und K- Halogenide mit F-Zentren scheint die Emission durch den inneren 
Prozess bedingt zu sein, wie die Messergebnisse zeigen. 

Um die Diskussion von Seidl besser vorfolgen zu können wurden in Fig. 3 die energe- 
tischen Schematen der Mechanisme wiedergegeben. 
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Photostimulierte Coelektronenemission 


Eine Möglichkeit der Haftstellenanalyse mit Hilfe der photostimulierten Coelektronen- 
emission wird an Hand der Abklingkurven während einer Beleuchtung mit Licht, sowie 
der Messungen von spektraler Quantenausbeute der zur Exoelektronenemission erregten 
Schichten von KCI zu diskuttieren versucht. 

Im Falle einer KCI-Schicht kann das auf der Fig. 12 angeführte vereinfachte Bandmodell 
betrachtet werden (Arsenjewa-Geil 1957). Dabei sind die folgenden Bezeichnungen zu be- 
riicksichtigen: 

ep — energetische Tiefe der F-Zentren, 

eg — energetische Tiefe der erregten Zustände, 

Ae — Breite des verbotenen Bandes, 

x — Austrittsarbeit vom Boden des Leitfähigkeitsbandes gemessen.! 

Wenn die F-Zentren in einem wirksamen Oberflichenvolumen (Oberflüchenschicht) 
von einer wirksamen Dicke generiert werden, so ist zu erwarten, dass ein Teil der Zentren 
nach der Beendigung der Generierung sich in einem erregten Zustand von der Tiefe ep 
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Fig. 4. Schema der energetischen Niveaus im Falle der verfãrbten Alkalihalogenidkristallen. V — obere Valenz- 
bankante, L — Boden des Leitfähigkeitsbandes. Über weiter Einzelheiten siehe Text. 


befindet. Dank der Warmeenergie des Kristallgitters (T~300°K) werden einige der erregten 
Zentren von den Elektronen entleert und die mit der Zeit abfallende Postelektronenemission 
in Dunkelheit beobachtet, wenn wihrend der Messung kein Nachfüllen der Chlor-Vakanzen 
mit Elektronen stattfindet. Wenn die thermische Energie des Gitters zur Emission der 
Elektronen beitragen sollte, so müsste der dem eiazelnen Elektron gelieferte Anteil etwa 
dem Wert der Austrittsarbeit x entsprechen, da die erregten Zustände sehr flach unter 
dem Boden des Leitungsbandes liegen. Der Wert von x wird für die Alkalihalogenide etwa 


1 Auch Elektronenaffinität bezeichnet 
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0,5 bis 0,8 eV angenommen (Wright 1957, Mott und Gurney 1950, Seitz 1940. Die Schitzung 
der thermischen Aktivationsenergie zur Emission auf Grund der Abhindigkeit (9) ergibt fiir 
T, =300°K * 

e~ 0,6 eV 


Dieses ist unter Voraussetzung geschãtzt worden, dass eine Glow-Kurve des KCI, wenn bei 
tiefen Temperaturen zur Emission erregt, ein Maximum in der Umgebung von Zimmer- 
temperatur aufweist. Die, auf diese Weise geschätzte, thermische Aktivationsenergie zur 
Emission von Exoelektronen bei Zimmertemperatur ist tatsächlich von der Grösse des 
angenommenen Wertes der Austrittsarbeit x. 

Die Kurve 4 in Fig. 5, die in doppelogarithmischer Auffassung dargelegt wurde, 
zeigt das Abfallen der Dunkelemission einer vorher durch Gasentladung erregten Post- 
elektronenemission aus einer KCl-Schicht, die an der Wand eines G. M.-Zählers nieder- 
geschlagen wurde (Sujak 1956). 

Die Kurve B in Fig. 5 zeigt die Abklingkurve der photostimulierten Coelektronen- 
emission, die nach dem Abklingen der Dunkelemission (Kurve 4) von der gleichen Schicht 
gemessen wurde. Das eingestrahlte Licht war von etwa 5700 A Wellenlänge. 


R GDG 


10 Cu lip D LN 
RN ID NR AUN LEN SIS EUN NS MAT MONTE 
Fig. 5. Abklingkurven der Exoelektronenemission von einer KCl-Schicht (eigene Messungen). 4 — Abkling- 
kurve der Postelektronenemission (im Dunkeln), B — Abklingkurve der photostimulierten Coelektronen- 


emission unter Beleuchtung mit Licht von etwa 5700 A Wellenlänge. Kurve B wurde aufgenommen nachdem 
die Aufnahme der Kurve A beendigt wurde. 


Aus dem Übereinstimmen der Gestalten der Abklingkurven A und B in Fig. 5 ist zu 
erwarlen, dass, wenn der Erregungszustand der F-Zentren (mit Elektronen gefüllten Elek- 
tronenhaftstellen) z. B. durch eine Beleuchtung mit Lichtquanten aus dem Bereich des Licht- 
absorptionsbandes der F-Zentren erzeugt wird, das Afballen der photostimulierten Coelek- 
tronenemission auf die gleiche Art wie die in Dunkeln gemessene Postelektronenemission 
verlaufen wird. Der Unterschied wird sich lediglich in den Intensitüten der Emission 
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zeigen. Dieses wird nur in dem Falle gelten, wenn die gleichen Haftstellenarten in beiden 
Fallen zur Auswirkung kommen. 

Wenn das eingestrahlte Licht Quanten enthãlt, die — wie das in Fig. 6 gedeutei wurde — 
nur ein diskretes Haftstellenniveau bedecken, so miisste eine Abklingkurve von der Art 
erwartet werden, die bei den Erwärgungen nach Nassenstein (1955) durch eine Gleichung 
von dem Typus (1) beschrieben wäre. Es wäre dabei vorauszusetzen, dass kein Nachfüllen 
mit Elektronen und Generieren der Haftstellen während der Messung stattfindet. Wenn 


"n Lichtintensität 


a) b) 


Fig. 6. Links: Energetisches Schema der mit Elektronen besetzten Haftstelle. Rechts: Intensititsverteilung 
des eingestrahlten, weissen Lichtes der Quantenenergie nach. Wirkungsvoll für die Emission können also nur 
die Lichtquanten sein, für welche h,> Ep ist. 


dagegen mehrere diskrete Haftstellenniveaus durch das eingestrahlte Licht bedeckt werden, 
so wird eine Abklingkurve von der Art einer Summe 


N 
yd Du Aj e- "it (19) 


zu erwarten, die sich; im Falle grósserer Unterschiede in der energetischen Tiefe der Haft- 
stellenniveaus sehr leicht bei der halblogarithmischen Darstellung In N/t, t auf die einzelnen 
Abklingkurven trennen lässt. 

Im Falle, wenn die Haftstellenniveauss nicht diskret liegen, aber eine Verteilung der 
energetischen Tiefe nach zeigen, wird vielmehr eine Abkliugkurve der photostimulierten 


Coelektronenemission von der Art 


= ER (20) 


zu erwarten sein. 
Auf diese Weise kann eine Haftstellenanalyse der Abklingkurven bei festgelegter 


Temperatur auch unter Beleuchtung durchgeführt werden. Dadurch werden gróssere 
Intensititen als im Falle der Postelektronenemission in Dunkelheit erziehlt. Wenn das 
eingestrahlte Licht tatsächlich nur einzelne Haftstellenniveaus bedeckt, wird sogar kein 
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monochromatisches Licht dazu gebraucht. Die Voraussetzung ist, dass es Quanten in ge- 
nügender Anzahl enthält, deren Energie Av gleich oder grósser als die optische Aktivierungs- 
energie für eine Überführung des F-Zentrum in den angeregten Zustand ist 


hv 2 £g €g (21) 


Da das Energieniveau des angeregten Zustandes angeblich 0,1 bis 0,2 eV unter dem Boden 
des Leitfãhigkeitsbandes liegt (siehe auch Bohun und Dolejsi 1959), kann folgendes ange- 
nommen werden (für Raumtemperaturen des Kristalles): 

Ept Eg~EF (22) 
Aus diesem Grunde kann bei vielen Untersuchungen weisses Licht benutzt werden, unter 
der Voraussetzung, dass energetisch erforderliche Quanten enthalten sind. 

Eine Analyse der Abklingkurven unter Beleuchtung gibt also auch die Möglichkeit, 
die energetische Struktur der Elektronenhaftstellen zu untersuchen. Die Ergebnisse einer 
solchen Analyse müssen jedoch sehr kritisch betrachtet werden und immer nur bei den- 
selben Versuchsverhältnissen untereinander verglichen werden. Die Änderungen der 
Austrittsarbeit für das Elektron sind nämlich sehr von den Zustandsänderungen der Ober- 
fläche abhängig (Wolkenstein 1960, Simon und Suhrmann 1958). Dabei ist auch zu berück- 
sichtigen, dass die Austrittsarbeit für die Halbleiter, je nach dem ob optisch oder thermisch 
ermittelt, folgend definiert wird: 


Poptisch — Ep +x (23) 
heka û ١٢٢ 9, 

@thermisch — 3 Ep PR (24) 

Wenn weiterhin für x positive wie auch negative Werte zugelassen werden, so kann die hier 

mit p bezeichnete Austrittsarbeit mit der vom Seidl eingeführten Bezeichnung e, identifi- 
ziert werden. 

Wenn die Wellenlänge des beleuchtenden Lichtes geändert und dabei die gemessene 

Elektronenintensität auf eine gleiche, auffallende Quantenzahl (oder Energie) reduziert 

wird, so ist im Falle einer selektiven Absorption des Lichtes K durch die mit Elektronen 


besetzen Elektronenhaftstellen (z. B. F-Zentren) eine dem entsprechende, selektive Kurve 
der photostimulierten Coelektronenemission zu erwarten: 


NT 


NK) (25) 


wobei E — Energie des in der Zeiteinheit auffallenden Lichtes, 
K — Absorptionskoeffizient, 
t — Zeitdauer der Messung bedeuten. 

Im Falle der Alkalihalogenide z. B. KCI oder NaCl, wurden unabhängig von Bohun 
(1960), (siehe auch Bohun und Mitarbeiter 1960) und mir (Sujak 1953, 1956, 1957) dies- 
bezügliche Messungen durchgeführt. Der Verlauf der Kurven für NaCl und KCI ist in Fig. 7 
und Fig. 8 zum Vergleich wiedergegeben werden. 

Da die Analysentechnik der Absorptionsbänder der F-Zentren für Alkalihaloide sehr 
gut ausgearbeitet ist, können mit Hilfe des Vergleichens dieser mit dem Verlauf der auf das 
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auffallende Licht bezogenen Quantenausbeute weitere Erkenntnisse über die Rolle der 
Oberfláchenzustande bei der Emission von. Exoelektronen gewonnen werden. 

Der Einfluss der Oberflüchenzustinde scheint sehr gross zu sein, da die genannten 

- Selektivitüten bei NaCl und KCI bis jetzt nur im Falle der Emission in einer Gasatmosphäre 

gemessen werden, (Sujak, 1953, 1956, 1957, Bohun und Mitarbeiter 1960, Bohun 1957, 

Asmus 1936). Im Hochvakuum sind diese Selektivitiiten nicht zu finden (Fleischmann 1933, 

Apker und Taft 1950, 1951, Klein, Schatz und Seeger 1960). Es ensteht also ein Problem 


und zugleich eine Möglichkeit, auf Grund dieser Unterschiede den Einfluss von Oberflüchen- 
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Fig. 7. Verlauf der optischen Absorption (Kurve a) und der photostimulierten Coelektronenemission (Kurve b) 
für KCI-Kristall. (nach Bohun und Mitarbeitern 1960). 
Fig. 8. Verlauf der photostimulierten Coelektronenemission für KCl-Schicht, die auf gleiche auffallende Licht- 
energie normiert wurde (eigene Messungen 1957). 


zustinden mit Hilfe der Exoelektronen-Methode niher zu studieren (Roubinek 1960). 
In diesem Falle muss aber auf die Untersuchungen der photoelektrischen Emission aus 
Nichtmetallen im allgemeinen hingewiesen werden (Arsenjewa-Geil 1957, Wright 1957, 
Mott und Gurney 1950, Seitz 1940, Joffe 1956). Ein grosser Teil der photostimulierten Exo- 
elektronen-Emissionsphänomena kann nämlich als photoelektrische Emissionserscheinungen 
aus den Halbleitern oder Nichtmetallen anerkannt werden. Der grundlegende Unterschied 
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in diesen Füllen liegt vorwiegend in den gemessenen Stromdichten und den Versuchs- 


verhãltnissen. 
An dieser Stelle muss noch darauf hingewiesen werden, dass die Potentiellschwelle 


(Energiebarriere) für Elektronen und die energetischen Bãnder in den Oberflächenschichten 
eines nichtmetallischen Kristalles nicht so einfach verlaufen, wie das in der ersten Näherung 
z. B. vom Seidl (1959) und von mir in diesem Abschnitt diskuttiert wurde. 

Es müssen weiterhin die Einwirkungen der Raumladung und der Oberflüchenauf- 
ladung, die während einigen Erregungsprozessen der Exoelektronenemission zustande 
kommen, auf die Neigung des Leitfühigkeitsbandes berücksichtigt werden. Dieses wird 
nümlich eine Móglichkeit des Entweichens eines Elektrons aus den inneren Schichten der 
wirksamen Oberflächenschicht anderen. Dieser Einfluss der Raum- oder einer Oberflächen- 
ladung muss sich besonders bei den Messungen der photostimulierten Coelektronenemission 
bemerkbar machen. 

Zu diesen Fragen wurden jedoch auf dem Gebiet der Exoelektronenemission bis jetzt 
wenige Versuchsergebnisse geschaffen, deswegen kann dieses Problem hier nur angedeutet 
werden (Matyás 1957; Hanle, Kanzler und Scharmann 1961). 

Für die photoelektrisch gemessene Austrittsarbeit kann, wie bekannt geschrieben wer- 
den: 

Poptisch = £p + % (26) 


wobei diese fiir das Elektron aus einer Haftstelle von der energetischen Tiefe e, in einem 
Halbleiter deren Oberfläche elektrisch neutral ist, gemeint wird (siehe Fig. 9a). We enn jedoch 
eine elektrische Überschussladungdichte auf der Oberflãche zugelassen wird, so werden 


গল - چم‎ ++ 7 
0 x Dat x een 
ù x 
pt = ae +6 u.a 5 
o p opt Taf dan +৫০ চো, ae 
a) bl oot" at Ade + Ê, 


c) 
Fig. 9. Biegung der Bänder als Nachwirkung einer Aufladung der Oberfläche: a) Oberfläche neutral, b) Ober- 
fläche positiv aufgeladen, c) Oberfläche negativ aufgeladen (Vergleiche Wolkenstein 1960). 


die Energiebänder (Leitfähigkeits und Valzenband) je nach der Ladungsdichte und Zeichen 
nach oben oder nach unter in der Oberflüchenschichtdicke von 1 gebogen. Dies ist in der 
Fig. 9b und 9c schematisch gezeigt. Die Biegung der Bänder ist also zu dem Fall nach 
Wolkenstein (1960) analog, wo die chemosorbierten Molekülen eine Überschussladungs- 
dichte an der Oberfläche des Halbleiters bewirken. 
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Wenn also das emittierte Elektron aus einer Tiefe J heraustritt und die Oberflãche 
eine positive Uberschussladung aufweist, so wird eine Herabsetzung der Austrittsarbeit % 
um einen Wert Ax erfolgen: 


@optisch = ê» + x — Ax (27) 


-Die Tiefe Z soll z. B. nach Wolkenstein (1960) um 10-4 bis 10-5 cm betragen. Dieser 
Wert entspricht merkwiirdigerweise z. B. der optimalen Schichtdicke der Exoelek- 
tronenemittierende Oxyddeckschicht, die etwa 107% cm betragen soll (Müller 1957). 
Es kann also angenommen werden, dass im Falle der Emission von Exoelektronen die Tiefe 
von 10-4 bis 10-5 cm die optimale Emissionstiefe ist, dass heisst die Elektronen werden aus 
dieser Tiefe am leichtesten ausgesandt. Wenn aber das der Fall ist, so kann erwartet werden, 
dass wenn die Oberfläche mit einer, der Grösse nach, gleichen aber negativen Überschuss- 
ladungsdichte belegt wird, erhöht sich die Austrittsarbeit um den gleichen Wert Ax für 
jenes Elektron aus der Emissionstiefe /. Die Austrittsarbeit wird in diesem Falle also be- 
tragen: 


Poptisch = Ep + 4 + Ax (28) 


Es wird auch verständlich, warum für die gleiche Lichtquelle und Oberflächengrösse 
öffters Photoelektronen aus den Staubkernen mit F-Zentren, die positive Oberflächen- 
ladung zeigen (nach der Theorie der Aufladung durch Brechen wird die Ladung an der 
Oberfliche lokaliesiert — Szaynok 1960), gemessen werden als aus denen, die negativ auf- 
geladen sind (Sujak und Szaynok 1961, Szaynok 1960, Sujak — Habilitationsschrift). Dies 
wird sich als gróssere Emissionsintensitát bei der photostimulierten Coelektronenemission 
auswirken, wenn die vorher durchgeführte Verfirbung der Ausgangskristalle eine positive : 
Überschussladung in der Pulver-Wolke bewirkt. Es ist nämlich zu erwarten, dass die Licht- 


quanten von der Energie hv, wobei 
er « hy « eg d- X (29) 


eine photostimulierte Coelektronenemission besonders von den Staubkernen die ober- 
flächlich positiv aufgeladen sind, doch bewirken werden. Wenn die positive Aufladung der 
Oberfläche eine Herabsetzung der effektiven Austrittsarbeit um einen Betrag 4x von der Grösse 
|Ax|>|x| bewirkt, so kann für die Elektronen aus der optimalen Emissionstiefe l 


folgendes geschrieben werden: 


@optisch < £5 (30) 
oder wenn für Poptisch die Bezeichnung ৪০০ gewühlt wird: 
৪০০ Sê» (31) 


Das bedeutet aber, dass die Exoelektronenemission lediglich durch den inneren Prozess — 
eine Uberfiihrung des Elektrons in das Leitfühigkeitsband — bedingt wird. In diesem 
Falle spielt also die Austrittsarbeit gar keine Rolle, wie das auch ófters beobachtet wurde 


(Bohun und Dolejsi 1959, Sujak 1956, 1957, 1959, Habilitationsschrift). 
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Aus der Ungleichheit (29) und aus der Fig. 9 kann auch verständlich sein warum für 
die gleiche Lichtquelle ein aus den mit -Strahlen verfirbten KCl-Kristallen hergestellte 
Pulver fast keine Photoemission zeigte, solange die Pulver-Wolke eine gróssere negative 
Uberschussladung aufweiste (Sujak und Szaynok 1961, Szaynok 1960, Szaynok — Disserta- 
tion). 

Vielleicht kónnte dasgleiche auch die Glow-Kurven zutreffen. Dabei miisste aber die 
Gleichung (24) beriicksichtigt werden. Eine Diskussion muss leider zur Zeit noch unter- 
lassen bleiben, da zu dieser Frage auf dem Gebiet der Exoelektronen bis jetzt kaum Ver- 
suchsergebnisse geschaffen wurden. 
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SPIN WAVE THEORY FOR CUBIC FERROMAGNETICS 
I. DYSON'S MODEL 
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The principal features of the Dyson spin wave theory of a cubic crystalline isotropic ferro- 
magnetic are presented. 

In the subsequent part of the present paper a new spin wave approach is developed, 
based on Matsubara's technique. 


1. Introduction 


In his excellent analysis of the cubic ferromagnetic Dyson (1956a, b) formulated the 
strict theory of spin waves. As the spin wave states are non-orthogonal, Dyson developed 
a special formalism of interacting spin waves and proposed a method of obtaining the thermo-, 
dynamical quantities within the framework of non-orthogonal states. There are two types of 
interactions in the Dyson model of a ferromagnetic. The one, termed kinematical interaction, 
is due to the fact that spin S (S — spin quantum number) of each atom can be located in 
space in 25 --1 different ways only, so that the quantum numbers of the oscillator operators 
co-ordinated to the spin operators cannot exceed 25. The other type of spin wave interaction, 
termed dynamical interaction, results from the assumption of a special form of 
basic states independently of the Hamiltonian. 

In the first part of the present paper Dyson's theory is outlined in general. The slight 
modifications introduced concern the generalization of the theory to any type of neighbour- 
hood of lattice points and the derivation of correspondence between the oscillator operators 
and the spin operators directly from matrix elements for the indefinite metric operator. 

As the kinematical interaction at low temperatures is small, as proved by Dyson (1956b), 
in the second part of the present paper we establish another form of the statistical function 
using the orthonormal set of basie states in the reciprocal lattice. We then apply the Matsu- 
bara (1955) representation. The ladder diagrams obtained thereafter enable us to compute 
any contribution to the magnetization. 

In the third part, we derive formula for the mean magnetization for three types of cubic 
lattice. The formula obtained is identical with that given by Dyson (1956b). As to the cor- 
rection to the magnetization proportional to 74, we have computed only the one term 


(983) 
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independent of S. Computation of the full coefficient requires consideration of higher order 
diagrams. 

Our main aim was to give a method for computing various thermodynamical quantities 
easier than the strict but difficult formalism proposed by Dyson. The present modification 


of Dyson's theory satisfies this requirement. 


2. The Hamiltonian of a cubic ferromagnetic 


We assume the ferromagnetic to consist of a single species of atoms with any number 
of magnetic electrons. The atoms are located at regularly distributed cubie lattice points. 
We put h=1. The vector to any translation lattice point is given by 


IN*— 11 


j= j'a, + j*a, + Pa; =a; j<| নিলা, 2 


where @,, 02, Q} are the principal vectors of the lattice and N — the number of lattice 
points. We apply the periodie boundary conditions. The reciprocal lattice is defined by 


A N* (4,b! zi Ab? T 440?) = N% Abi, Ài < 2 , 2.2) 
i 1 ikl 


where & is the permutation tensor. 


To any lattice point j is attached the spin operator S; satisfying the commutation rules 


[SF Si] = iOS, [SF Se] = iê ıS, LS, SA] = 2৯০, (2.4) 
ep > uet: 
55 = ৪32, (2.5) 
whence 
نوم يدد :لس وول ال لل‎ Solis ê 
[5১] = 2৮১৮ 1555] ২75857৮1355] = 2989, (2.6) 


We write the Hamiltonian of the system as follows: 


8# - L 73-1 29(Jj-k 
J J#k 


er > $43 9(0) >} Sas +>,9( j-kp S; S,. 
j j jk 


mH 
5 


Le A (2.8) 
wherein S is the spin quantum number of each atom, H — the magnetic field strength, 
(|j kl) — the exchange integral and m — the magnetic moment of each spin. 

Since 


S? = S(S+1), ^: „dA 
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we have 


H=49(0) NS(S+1)+L ZS -4 NIG k S S. (2.10) 
J k i 


J» 


We can now expand the exchange integral in a series over the reciprocal lattice, 


; 1 B 
Ilj- k) - ( JU, 9ı - 9_,, (2.11) 
A 
whence 
1 
900) = 4 P3 I (2.12) 
A 
and, conversely, 
1 4 
== 2 (j — kj) e" a-». (2.13) 
J 


In the special case of interaction only between nearest neighbours we denote by 6 the 
vector connecting the nearest lattice points and get 


k=j+ó, 909৯9, 901)-7) = 9ó(k — j — à), 
1 A 8 zo PAS ۵ 3 এরি, ê ۹ 
gam (= 3 « 2 ১? 91 iiU- (lg — j - ê) = 
j,k Jk 
- 2 )) ود‎ y=) oki, (2.14) 


726 6 


The spin operator attached to a Ã point in the reciprocal lattice can be given by the 
Fourier series 


S,=N"h 2 e AS. (2.15) 
J 
By the last equation we can rewrite Eq. (2.10) in the form 


LEM Sie t] Dio E LNA 1819, S,- 6-5, (2.16) 
A A 


where we have used the Kroenecker delta 


600) = N3 2 dba, | (2.17) 


3. Spin wave states 


Let us denote the ground state of the system as |0> and define it by the equation 


S? [0»— — S0». (3.1) 
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From the relation 
S S aa + iê) - (ê - ıê) = ptr — UŚ, S = 
= $} - Gps  9)ق‎ +1) — 6n 3. (3.2) 
we conclude that 
SFS-|0» = S(S+1)[0>—(57*|0> + $|0>=0 


or 


Ś; |0> — 0. (3.3) 


S; is the spin deviation annihilation operator and Sj — the spin deviation creation operator. 
Let the complete set of orthogonal but, not normal state vectors be determined as 


|u> = u [25)- (uj) (57)5]]0 >, (3.4) 


wherein u; is the positive integer and [J denotes the product over all points of the lattice. 
J 

We now ER on the norm of the vector of Eq. (3.4). For this purpose we notice that 

S$; ($+) = (SPS — 2($7)8- 1[ pS? (51 pH1yp2, 3; 1 (3.5) 


The above relation can be proved by total induction on the ground of Eqs. (2.6). Eqs. (3.1), 
(3.3) and (3.5) now yield 


S (Sf) |0) = AS- | — Su; + (“| I» 


- Uj الس‎ é 
= 5۰ Uj (ı — “35 ) (5)4-519 Xê 


ক u; — 1 vam DN A 
(Š3)? (5910) = (25)? u; (u; — 1) Ê 25 )(ı ES ZZ 


and, quite generally, 


LOCS) "i (St) VOD 0 u; u;! F(u,), (3.6) 
with 
1 2 u— l 
We have, finally, 
< v|u = Ji 69০, Ti = II F(u,), BUSA ১০) (3.8) 


We now introduce the RON set of orthonormal SSE” state vectors 
18) 11101) 70/05]10), (3.9) 
j 
(এ (01171105071, (8.10) 


J 
[nj mal = [mj ايا‎ [nj mı] = 25. — (811) 
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z 


Let F be the indefinite metric operator in the space of states (3.9). We require that 


(v|F|u) =< v|u y (3.12) 
or 
(oF (৮1৮৮0117170) (v 7^], (3.13) 
J 
1 x TA 
Fo) = [I 1 = ag iin — n|: (3.14) 
=1 
We shall prove that 
U, TS 1 * Ki 
(01110) = (0) | (1 হত mm] | - (3.15) 
We have 
8 1 1 
|». 1- gg (nim — J = — 25 n» "nu = — zę > 
10. Im: 
n; |17-2507//-0 | = 11755256777 l+ 1) |n, 
vor TM 
„ji — gg (06% — J = | নু nin t 
yielding 


vj 
uf v 1 ts Û > j 
1170) = 11 E — agn J a |: E J i 


DN s Ipe Lat. A 
x [ı- g مې نه - وم‎ || as (Mi Ni »|-n- zema qug X 
1 * 1 * 1 * 
x |1—5sGi 71) |--|1— হুক 1-0 | =D 35 ail 
1 * US 1 * D 1 * 
x = nin | nj ` nj! des || nun] = 


il i v,—1 1 çê il " vj 
= | [1- چي‎ |n mI gD =|- 5s J1i Î - 


Thus Eq. (3.12) assumes the form 


017 le uj!) * I = = in n | “a, Û 


KOOS Fr (o1 wd” (5) FTA: (3.16) 
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Thereby we have established the following correspondence, Morkowski (1959), 


^ ; 1 
ر‎ > ( l- zg" „) UE (3.172) 
St > (2 Sjn, (3.17b) 
S — S4 N; Nj (3.17e) 


The last equation (3.17c) results from Eq. (2.6), 


In order to establish one to one correspondance between the space of physical states |u» 
and the space of ideal states |u), we have still to add the condition 


0 <u; <25. (3.18) 


The Hamiltonian (2.10) in terms of the ideal operators 17, n; is now 


S, S,- 1(6 مغ‎ Sz) کر‎ (3.19) 
[aee E,+ LX) njn;+ 1S ১9011 — kl) (e; — m) (n; — m4) + 
J Jk 
+} >; JJ - El) nj v (n; — M)? (3.20) 
j,k 
E, =— LSN+ 49 (0) NS(S--1) — 4 5? 51900) — kp. (3.21) 
Lk 


Taking into account the nearest neighbours only, we have by Eqs. (2.13) and (2.14) 


H = 974 2) nj Nj +495 >) (nj = 1:49) (0; — 040) 29 2 ni nj cl; — "49^, (3:22) 
y قرز‎ 7১9 
E, = — LSN — 9158০, (3.23) 


as 
IO = X3*19,- 9N- Y y,- 9N اک‎ Şj FFF y (6) =0, 
Ã A 0 A 0 


where yọ is the number of nearest neighbours. 


Obviously, we might also consider more highly complicated cases, e.g. including the 
exchange interactions between the nearest, next nearest, third, fourth and further neigbours, 
by putting 


(lj — Kl) = 9.8 k — j — 6) + 9,6 (k= —3,) + 9,5[(e —5 —à)--.., — (8:25) 


5 
0, are the vectors from a given lattice point reaching to its nearest, next nearest etc. neigh- 


bours. 97, i=2, 3, ... can be treated as variational parameters; they can be deduced from 
the minimum condition for free energy. 
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We now introduce the Fourier transforms (ideal spin waves) of the ideal spin operators 


gam 21 ent dry, (3.25) 
j 
=Nh DM qs (3.26) 
j 
with the commutation rules 
= [ar a - 0, .[e, c, = Oy n. (3.27) 


The Hamiltonian (3.20) transforms to the form 


Jk = Eo (L+ eoe — 16 ez 5 ai xo WALT (3.28) 
Ey =— LSN+ 1 S(S+ 1) > O, 9051৬, (3.29) 
7 
و‎ > a! (I= টা, dict Dts) (3.30) 
= (Fo - 99). (3.31) 


The non-diagonal put of the Hamiltonian (3.28) describes the dynamical interaction between 
the spin waves. JI, are the coupling functions of spin wave interaction. 


4. The statistical function 


In full accordance with Eqs. (3.17a—c) and (3.18) we define the statistical function as 


Z = Sp exp (— BH) == 222 (u| exp (— BL) |u), (4.1) 


wherein 
1 MEROE < 25, 
p a T E, = IT Eu; Eu, = 8 u; "sd (4.2) 
J 


In formula (4.1) summation over u ranges over all numbers u, of the state vector |u). 
Let us now introduce the set of orthonormal ideal spin wave state vectors 


lo) = Ta) «510: (4.3) 


a, is a positive integer (zero included). The product [runs over the entire reciprocal lattice. 
a 


By Eq. (4.3), we have 
Z চা E, SJ (ula)(a|e""7|b)blu). (4.4) 


a,b 


Non-vanishing matrix elements in the above expression occur for 


a = AB 22/7 o GEM, ZonN. (4.5) 
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We split the Hamiltonian (3.28) into three parts, 


2 سن‎ Ry (4.6a) 
Hy = Y1 (L+ eoe; (4.6b) 
A 
Hy ko 19 N pS bevel Meg NC OT (4.6c) 
290 


With respect to Eqs. (4.6b, c), we have 
e- ه هو‎ — e-B%S(8), SB) = eB He eK Pot PD, (4.7) 


The operator SE) is given by, (Matsubara 1955), 


ê B 8 
909) = e = | 0 | 07৮7 [9105) ... 210০], (4.8) 
n=0 0 0 


n 


wherein the dynamical operator 


H(t) = er Hope "Fe (4.9) 
is written in Matsubara's representation, and P is Dyson's ordering symbol 


71%105)%1(0)...910০)] = H(t) r(t;)..<#r(T;,), 


TA 72725 ST 
fı و1‎ 


(4.10) 


in’ 
m A T. , . ^ . . . 
[he operator P may be replaced by Wicks's ordering symbol 7, since both are identical for 
the boson field. 

The statistical function (4.4) is now 


Zoe PPS E SJ (ula)(ale-*%*Š(8)|b)(blu). (4.11) 
u a,b 
The state vector (a| may be written in Matsubara’s representation 
Ua» = (a (B) |S (8)! b) . (4.12) 
wherein 
(a(B)| (01111702)(%)1) *]. (4.13a) 
Ã 
| = II [(5, 9) -(25)^4]|0) (4. 13b) 
and 
&,(T) = e” Hoa, en Ho = e" 2 ০), (4.13c) 
a(t) = e Figi ec ns = ০৭০) 0, (4.13d) 


We can now apply Wick's theorem to the expression aj (ti)æ (T2). We have 


fa; teast] = N [25 (5) a (JI Ta (r) a). (4.14) 
E 2 
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We shall term the last term in eq. (4.14) „z-contraction”. N is the symbol of the normal 


product. The arrangement of operators in the N-product is the following: all annihilation 
operators stand to the right of the creation operators, i.e. 


Nase) (ta) = a(t) (3). (4.15) 


The T-product has the form 


may | CD (rg), T, > Te, 
T و و‎ 2 : 
CACAI | ae (4.16) 
whence 
* E > Tı > T9 > 
OC রো তীর 
Since, according to Eqs, (4.13c) and (4.13d), we have 
a.) en (2))] = ne îron [ayê] Oge DD, (4.17) 
the t-contraction now becomes 
* | 0, T1 = T2, 
& (T4) & (5৪) = | বা (4.18 
1 1) ( 2 UC oda aX 1 2) Ti ES Tos ) 
or 
o; (5) &,(T,) = Ö, o B (ts = Tı) e ran, (4.19) 
| | 
wherein 
٢ TU, (4.20) 
Ple) i= l0, t= 0, 
0 = z | da >, Del: (4.21) 
2771 % — 16€ 
Quite similarly, we obtain 
%a (Ty) a (Tą) = Ö, oð (v — Ty) "rao پا‎ গা (4.22) 


We can now provide a graphical interpretation of the consecutive terms in the S(f) 
expansion. With respect to Eqs. (4.8) and (4.12), we have 
B 


B B 
Uas = (a(8)|b) — jl ৫5০91910919) + >, i dr, ji dz, (4(9| T [Fr লে) x 
0 0 0 


B B B 
x x (rJ) — 3 Í dr, f dz, 1 de, (a (91771002150) +... (4.23) 
0 0 0 


The first term in the foregoing series represents free motion of the q= 2) a, spin waves, 
7 


8 
(a (6)| b) = PODR : Öz LLE 0 L 
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The second term in (4.23) can be written as follows: 


(a (8), (21b) --15N-' DTA, (Ql JT (as! 6,17% aż (B) ৮০4১0) x 
200 = 


x a 4(T) a,(®) ০০(0(৮]10), (4.24) 
The non-vanishing matrix elements satisfy the condition 
DLLODLPP 2 (4.25) 
p n 
and the Feynman diagram is of the form: 
d 
U, |, [us و ولا‎ 
0 


Fig. 1. Diagram illustrating the matrix element 


(a(ß)|as+ a(t) ag a(t) aq(7) qu(z)|b) 


We now apply the Wick theorem to the third term of (4.23), obtaining a single N-product 
with all operators, eight terms with N-products times one T-contraction, and four terms 


with N-products times two T-contractions. The normal product with all operaiors yields 
two disconnected diagrams of the first order, as shown in Fig. 1, whereas e.g. a term with 
one T-contraction yields the following graph: 


g 


0 
Fig. 2. Diagram illustrating the matrix element 


do 1000) ag, (r3) (a(B) |g, — 1, (v1) oo, 1, (0) 09, — 2. (79) Gg Tya a, (T,) 2a, (72) |b) 
re ê سنا‎ || 


8 


0 
Fig. 3. Diagram illustrating the matrix element 


+ 102) (Ta) dg, — a(t) “a, (Ta) (a(8)| a, a (73) Que 8,022) ap (1) ag,(T,)|b) 
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The Feynman graph for two t-contractions and the corresponding matrix element is 
represented in Fig. 3. 


The technique of constructing an arbitrary diaeram is resented in Dyson’s paper 
q g y ৪ p y pap 


(19562). 
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SPIN WAVE THEORY FOR CUBIC FERROMAGNETICS 
II. NEW APPROACH 


By JAN SZANIECKI (ŚLEDZIK) 
Institute of Physics of the Polish Academy of Sciences, Department of Ferromagnetics in Poznań. 
(Received May 31, 1961) 


A new spin wave-approach to the theory of a cubic ferromagnetic is developed. The 
smallness of the spin wave kinematical interaction at low temperatures is taken into account, 
which enables us to compute the statistical function as the trace over the complete orthonormal 
set of ideal spin wave vectors. The Hamiltonian of the ferromagnetic is taken from  Dyson's 
theory, but, owing to the new way of obtaining the statistical function and applying Matsubara's 
technique, the logarithm of the grand partition function is now virial expansion containing 
Bloch's term and the sum of connected ladder diagrams. 


1. New definition of the statistical function 


In the first part of the present paper we referred to Dyson's definition of the statistical 
inction for a ferromagnetic at low temperatures. In his second paper Dyson (1956b) estimated 
ıe kinematical spin wave interaction to be proportional to exp (-aT,/T), where a(>0) 
a number of order unity, 7T — the absolute temperature, T, — the Curie temperature. 
1 this way, kinematical interaction obviously being a purely statistical effect is quite small 
| suitably low temperatures. This fact enables us to proceed to a new definition of the 
atistical function which, following Bloch (1930), we define as 


25 
Z = Sp e" 5% =); ১ 0 (১10) — q) (ale $7 |a), (1.1) 
4=0 a 7 
here, according to Eq. (I 3.28) (part I of the present paper will be henceforth referred 


) as I), 
ż= Eg DA (47 a EINE le وه‎ pa ০2০০ (1.2) 
7 80০ 


yson's definition of the statistical function 
Z= DE, (u|e-?* lu) (14.1) 


troduces certain restrictions as to the localization of the ideal spin waves around a reciprocal 
ttice point. This becomes evident when we replace the oscillator operators n,, 7; in |u) Eq. 
(995) 
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(14 4.1) by their Fourier transforms &;, &; by means of Eqs. (I 3.25—26). The sums of factors 


exp 01 j,) appearing thereupon restrict the number of spin waves that can be localized in any i 
point. Since, as stated above; the kinematical interaction is small, we may neglect these 


restrictions and assume q<2SN spin waves can propagate simultaneously through any À 
point. It will be shown in Part III of the present paper that the new definition (1.1) of the 


statistical function leads to the same results as obtained by Dyson. 


As in Part I, we expand the operator exp (-PX), 


ZEE (1.3)‏ و 
١٢ (1.4a)‏ و9 2 )1 EF, =— LSN+$ S(S+‏ 
Hy = 2 (L+ eta: (1.4b)‏ 
Hy = — 19 N31 31 P5 uo an amas (1.40)‏ 
200 
in the perturbation series‏ 
rê 8‏ 
e-b% — e- 8E oP S (8), S(B) = DL 2l dry... | dz, Tir (0)...%1 (<), (1.5)‏ 
0 0 
with‏ 
)1.6( مو H(t) =e UH‏ 
in Matsubara’s representation.‏ 
Equation (1.1) now assumes the form‏ 
2SN |‏ 
gla|e""FS(B)|a). (1.7)‏ — په Z‏ 21402 :£56,752 
q a‏ 


On expressing S(8) by its perturbation series and on retaining the first term only, we obtain 


Bloch's partition function 


2SN 


Z9 = ৮৮ 3 Y ê (Ş) a, — q) (ale? 749) 
Ã 


4=0 a 
We have further 


٤۴) ٢٢ !)-*] =‏ سو ٢ل‏ پر 


= exp [- f 2 (L + s,)a,](0| 11 [xa (a) 2] = exp [= f 2 یا‎ ez)a,] (a), 


whence 
2SN 2SN 2SN 


£o eh es 2 ره‎ a 


9=0 a,=0 wm. 


(1.8) 


(1.9a) 
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The difference between Eq. (1.9b) and the expression 
(oy E SA R 
مه‎ উঠার PRZEJ PE لوي سه‎ (1.10) 
45,—0 Ay 4 
is of the form 
, iia 
200) _ 2700) — 2—8E ۳ e 48L = e abe | > سف شا‎ ten an 
q=2SN+1 7 ورت‎ 


+ e Bler (Dea 4 e-Pla-Deyt2er) 4 سی‎ Bes ta Der) 4 


+ e Ma Beuta) + e> Bear 73e... ] + > [০-%৫4-25741541)4- 
hı< kj Ak 
+ ePlemta-Deyter) 1 99667475747 ৫-2)5%) 4 e Plla—3)84,262,+824,) 1 


+ 8-8৫৫-3)5/4757/7857%) |. ړه8)2 سی‎ t+ —3)e gy 634) 4-6-96/44-4-3)৪/4+ 859) ai 
$e Partt وول رد پس‎ E S +...) (1.11) 


Even at zero magnetic field strength, the condition 8>>1 (low temperatures) and the presence 
of at least (2SN+l)e,,, (A44 <... <A), terms in exponents guarantee that the difference 
is negligible. 

We can now write 


cM > (a|e 9^ S(B)|a) 
Zo = <S(P) >= > (ale-#P.Ja) 


(1.12) 


Z = ePE y (a| e-?$*| a). (1.13) 


a 


Summation over a ranges over all positive integers a;— 0, 1, 2, ..., + oo. 


2. Matsubara's representation 


In order to compute the average < S(f) ১, we introduce the formalism developed by 
Matsubara (1955). We have 

a,(t) = ehe Po = ge T (2.1a) 
a ir) = et Hoge tHo = ge tt), (2.1b) 

We now split the operators %,(?), «(7) into positive and negative parts, 
a(t) = a (e) Fa), a(t) = eS) + a), (2.2) 
28১05) = وو‎ Ter. CZU) 5 78 (2.3) 

0 = (1 - gda EHD, اړو‎ = (17980 


The coefficients q} and w; will be determined later on. As usual, we define the normal product 
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| 
of a certain number of field operators as a product in which all negative operator parts stand 


to the right of the positive ones, e.g. 
Nja T) ০2] = N flag" (r) + af (zı) [as م)‎ + as?a] = 
a Nie (P(r) a; a (7) Pad (zg) ać (a) ړب‎ (x4) CA toj (ty) a ?(12)] = 
= a (x) 25 (ra) kal (xy) a (Ta) + az C (vg) a (v) + af (14) a(t) = 
= (14) 4(7) — [af (rı, a (T2)] 2052) (tq) 41205), a (13)] — 
om [aS (z,), ১4811) (3.4) 


The Î-product of C(7 (G(T) is 


1 دد‎ = Gato. سي‎ (2.5) 
The r-contraction is now 
fase) ar] = Nl) aid] + a(t) arg. (2.6) 
a (n) (ta) = 2 e হাহ Kara DN ٢ TEN eu 
or, with. respect to IL Ya (2.1a, b) and (2.3) 
„-(L-+ej)(zı- 2 
2H = 9 ay Tere) A EU 
We now require that the average of the N-product shall vanish: 
î Zi (a |e- ^N [az (tı) a$ (Ta)] |a) 
Mate) ab (r) = سي‎ ES 
This yields : 
€ %)(5)%০(2) Ex [a (z,), ০০৫+)(2)] ) = e Eten . g(L+eo)ta ১৫ 
x[< ao; » — 0,,591(1 = 0,)] = RT N ahia w 
x [< eı > — q,(1— ০9] = 
(am ) =< ata) +l = m+ l = qil — 0), (2.9) 
n, [69৫৮৯ 1]-1 (2.10) 
and 
afty) (ta) = 0,90 EED 09) — r)(n. + 1) + O(a - zı) nıl, 10 
> Û 


Similar computation yields 


পা) sr) = 6০০০৫০৮1905 — zı) nı + O(a (ng D] 218 
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and 
— w,(1 — q,) = m. (2.14) 
By Eqs. (2.9) and (2.14), we have 
77271 lc e[n;(n, +1)]%, 0,= — n; + e[n;(n, tl], e=+1. (2.15) 


In the special case of coinciding variables r, we obtain 


a,(t)@,(t) = 0, ,(n,--1), (2.162) a;()m,(r) = ûj, a. (2.166) 
[ | | | 


Direct computation also yields 


< N[aã(zı) ০505) %(T3) % „(va)| > = 0 (2.17) 
and, quite generally, 
< Ne; EDS) es NEA CZA (T+ 1) Zo (Tn+2) e Uu U 
np IM PESO, (2.18) 
as was proved by Thouless (1957). Thus, the T-product may be written directly as the sum 


of products of suitable number of T-contractions, which, being c-numbers, can be put outside 


the symbol of averaging. 


3. Graphical interpretation 


By (1.12), we have 


+00 0 B 
AOT ৯,৮2৫), = | a... | ctf a ite 
n=0 5 . 


0 
B 


19 B 
=Z0 fh 1 dz (((و)ر#) و‎ + xf is | dr, (T [Xr 0) ZN] > — 
۷ 0 0 


8 8 B 
E E dr, Ji dr, 1 dv, (T لج)ر#)‎ r (2) r (13)] > hak (3.1) 
E 0 0 
with 
H(t) -—19N93 XV iuo; (o) as). (3.2) 
hoo 


The first term in the expansion (3.1) now assumes the form 


B A 0 دم‎ * * 
rı-- | de< [ŻA] -39N3 DMA, | de CP [ae ae) )په 2 په‎ 2.(2)]> = 
5 Aoc 0 


-$ 9141 ১715 if dr [€ a, a(®) %,- 4(7) (x) oy > + ৫5+7(5) %-705) %ę(7) ০০()৯] = 
E Ê PSE | | | 
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دس وس TÖz‏ نس و , , سی Ze fies anka‏ 19-1 


= t 99 D NES ET +4 4 SWP p? I oqo‘ (3.3) 
20 0,0 
The Feynman graphs for both cases are illustrated below: 
o+ 9-9 esa 
o o 90 
aê alt) a5 — (<) aglr) ৫০() a. (7) aş. (1) )یه‎ a(t) 
! ! | | ! 1 | I 
1৮৮৬1 1 | 
pie 


Let us now analyse the second order term in eq. (3.1): 


1 A W Ice Nê l 
2T fax fan TlH (rı) #ı(r»)] = oT FZ Pa heê X 
Zm m 
B B 
X Î « [aux Tose, (Tı) Ge (Ty) Tę (Ty) a (T1) ০৮৮+৮ (Ta) Ku. (T2) Xu (Ta) ©, (Tall ১. 
0 0 


t-contractions between the same r-variables yields the disconnected diagram of the first 
order equal to 1/2! f. We consider therefore the graph with contractions between different 7. 


We have 
1 1 ê Ã x 
n= (ion ) Nue 
6و‎ 
xuv 
8 B 
x [dr [de [i a(t) a, (rı) a,(rı) a,(rı) aera) জপ পে এ+ 
0 0 — 1 ! | 


eset] | | | 
L3 1 | 
| 


* 


3 ৫৫০ (©) 2৫70) : چا‎ % (T1) a, (Ta) a p= „(T3) & (T3) a, (T2) + 
i | | | 
| | | 
i | | 


| l 
| 


p c, zu "E (n) a (11) % (r1) "x. „(Ta) NAM % „(Ta) %,(T») + 
4 | | | | | | 


E s cu EE i | 


! 
| a OWĄ 


= 


* M (n) (ai) at (T1) a (1) EAS CN PMI & „(Ta) থ,(22)], 
| 
| 
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By Eqs. (2.11) and (2.13) we get 


8 
lc 
Is = 21 (a gN-1 ) 2 103 re, | | ar. elut ty 9 — 8g) (T1— Ta X 
= a 
x [®(zı- Ta) N „n,(n, 351) (n, +1) m A(T, —T) (n. ma) (n, rl) nons] s [99 TAE »)ó(0 TE Ne M) x 
x ê(9-+ x 5) 0(u— x—0) + (65-42-40) 6(0— A— 9) I + 2-0) 0 — 5-2) + 6(6+4— 5) x 


x 6(6—A— n) 8(9-- ou x— 6) +0(6+A— h)6(6-1—5)6G6+%— 8) dd — x — 9. 
(3.4a) 


Integration over 7, and Tą yields 

6 8 
277৮৮৮৮৮৮21 
0 ৪ 


8 Ti B Ta 
Ź Tha] = far eat | dr, e” nn, (ng + 1) (n5 + 1) + IE [4 en x 
9 9 9 0 


x (na + 1) (n, +1) nons = £ [(" +1) (n, +1) ngno — nun,(nę +1) (na +1)] + -ٌ [(e® 1x 
x nyn,(ng + 1) (ny +1) + (e-*? — 1) (n,+1)(n,+1)rnene]j, a= ৪৮7৪, Sg — یه‎  (3.4b) 


Albeit, from Eq. (2.9) we know that 


De MEI مو‎ (m+1) („+ م46 لل‎ 
e ——, efa = : 
ną ny My (ng 4-1) (no 4-1) 


ną = [89৮৯ — 1]-1, يو‎ = 


whence 
(e® — 1) nn, (n, + 1) (n, + 1) + (e 9 —1) (n, + 1 (n, +1) nan, = (n, 4-1)(n, + 1) nin 
— nyn,(n, + 1) (ng + 1) + n,n,(n, + 1) (2, F 1) — (n, - 1) (n, + 1) nn, =0. — (3.40) 


Making use of the Eqs. (3.4a, b, c), doing some interchangements with vector indexes h.v 
and taking advantage of the symmetry property 12, = 157, 12, = T° ° , we get finally, 


z ON? Ê ` TE : ; (Gig ED (ny + 1) x 


Eg ^l وه‎ Er Bp 


Auv 


5৫7৮4717777 > TATU (Na + 1) (nya + 1)]. (3.5) 
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One part of the ladder diagram of Eq. (3.4a) has the graphical shape represented below: 


a 


, ৮ 


0.3 دم‎ 
Ti 


do (rı) ap — lt) په‎ (m) ao (rı) y+ se (ta) rd x(ra) a, (a) a, (75) 
l 


| | | i 


| 


* * 
Got a (0) Qo, — 2, (71) a, (1) Gs, (14) a, 7502) 2 Gare J. (75) A, 72) ৫০ (To) ay, + + 7১03) on په‎ 13) ay, (Ts) Q, (T3) 
| l 
| | 
| | 
| BE | 


( 
| | | ne ai = | | 
| 


Fig. 3* 


We denote the connected diagram, which in general consists of several parts, by مل‎ 
wherin i points to the order of the diagram. We include in 7; only those part of the connected 
diagram that are distinct as to their physical meaning; thus we obtain 


^ 


+00 
4 S(b) ১ = eli that yt... 22 exp (> HY (3.6) 


The grand partition function is now 


1-00 
Z = 252 (D573): (3.2 


i=1 


* In Fig. 3 t, should be in place of Ta and inversely. 
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By Eqs. (14.3), (1.4b) and (1.13), we have 


20) = e" BE x (aje. Poja) = 6689 m3 (aa!) (0| IT [ate x 


a;=0 u 


x e Pa „0 1|0) = e BB J] ` ( a! ED HI [a] (B ) (e )“]|0) = 


TERE 
— e- PB J] د دوا‎ ৮57118727৮5 (3.8) 
à 2830 P 
We used the well-known relations 
a0), 1=0, (0e: = (0|. (3.9) 


Thus, the function Z finally assumes the form 


7-০০ 


Z =e "*J]|l=e PW] exp (2119 (3.10) 


À ZA 
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EFFECT OF pH ON THE FLUORESCENCE OF FLUORESCEIN 
SOLUTIONS 


By M. RozwADOWSKI 
Physics Department, N. Copernicus University, Toruń 
(Received June 14, 1961) 


The absorption and emission spectra, mean decay time, relative yield and anisotropy of ` 
fluorescence emission of a fluorescein solution of the constant concentration of 10-5 g/cm? were 
investigated as to their dependence on the pH of the medium. The absorption and emission spectra 
were found to vary according to the pH value. This effect of pH can be explained by the molecule 
in the excited state giving off a proton. 


Introduction 


Owing to its intense fluorescence in the visible region, fluorescein is one of the dyes | 
whose photoluminescence has been most often investigated. However, the fact that different 
ions of the dye, presenting different optical properties, are present simultaneously in the 
solutions throughout a wide range of pH values has not been taken into account in all cases. 
This is of especial importance in the investigation of the photoluminescence as dependent 
on the dye concentration, the temperature and the solvent. As these parameters are changed, 
the relative concentrations of the different ionic forms of the dye molecule can vary. In 
strongly acid solutions, the green-blue fluorescence of fluorescein is ascribed to positive 
'jons, whereas the green-yellow fluorescence exhibited in basic solutions is considered to be 
due to bi-negative ions (Batscha 1926, Levshin 1951, Pringsheim 1949, Fórster 1951). 
In certain neutral solutions fluorescein fails to exhibit fluorescence (Howe 1916). Zanker and 
Peter (1958) carried out a systematic investigation of the effect of the hydrogen ion concen- 
tration on the absorption spectrum. In order to obtain solutions of different pH, they added 
water, ammonia, or sulfuric acid to solutions of fluorescein in dioxane. They considered the 
spectrum obtained from a solution in a mixture of dioxane and water to be due to univalent 
negative ions, that of a solution in dioxane and ammonia to bivalent negative ions, and that 
of a solution in dioxane and sulfuric acid to univalent positive ions. Fluorescein in pure dioxane 
possesses no absorption band in the visible range; this points to the presence of the colour- 
less, neutral form in the solution (lacton form). Grigorian et al. (1960) obtained the emission 
spectra of similarly prepared solutions. By adding 0.1 n sulfuric acid to dioxane, they observed 


( 1005) 
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the appearance of a maximum shifted towards longer waves with respect to the band corres- 
ponding to bivalent negative ions, whereas when 12 n nitric acid was added the maximum 
appeared on the short-wave side. The authors considered these differences to be due to 
the effect of sulfuric or nitric acid. Like Zanker and Peter, they correlated the different 
maxima in the emission spectrum to the respective ions of the dye. 

As the solutions used by the foregoing authors in their investigations of the absorption 
and emission spectra were prepared in different conditions, the necessity was felt of a more 
detailed investigation of the effect of pH on the fluorescence of fluorescein by measuring 


the different parameters for one and the same solution. 


Preparation of the solutions, and methods of measurement applied 


Chromatographically pure fluorescein of the red variety was used. All solutions were 
prepared as follows: fluorescein was dissolved in twice distilled water with an admixture of 
ammonia (required for full dissolving of the dye); then in order to produce the required pH, 
NH, or H,SO,4 or HNO, was added. Nitric and sulfuric acids were chosen as the ions SO, and 
NO, are very weak quenchers of fluorescence (Pringsheim 1949). The effect of reabsorption 
and secondary fluorescence did not exceed the accuracy of the measurements. For all solu- 
tions investigated, the dye concentration was 1075 g/em?. The solutions were stored in dark. 

The pH values were measured with a battery pH-meter made by “Ridan”, Warsaw. 
A glass electrode was used. The measurements accuracy was about 0.05 units of the pH scale. 
As the hydrogen ion concentration changes with storage (especially in weakly acid and 
weakly basic solutions), the pH values were checked before each measurement. Fig. 1 
shows the absorption spectra of the same solution measured after different periods of 


storage. 


0.3 
0.2 


0.1 


390 400 410 420 430 440 450 460 470 480 490 500 310 

A (mu) 
„Fig. 1. Absorption spectra from one and the same solution, measured after different times of storage: after 
20 minutes (pH = 5.10) — — — — after an hour (pH = 9.42) —.— .— after 12 hours (pH = 5.77) 
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The absorption and emission spectra were measured by the photoelectric method, 
using a Zeiss mirror monochromator and EMI 6294. photomultiplier in connection 
with a galvanometer (“Skalengalvanometer”” Zeiss-Jena). When the temperature dependence 
of the absorption was measured, the standard and the one investigated solutions were 
placed in a thermostatized vessel yielding stability up to 0.1% C. 

The measurements of the mean decay times of luminescence were carried out by means 
of a fluorometer built together with R. Bauer (R. Bauer and M. Rozwadowski 1959). The accu- 
racy of the measurements depended on the intensity of fluorescence and exceeded 10-19 sec. 
In this fluorometer, an ultrasonic birefringence modulator of higher luminosity and greater 
simplicity than the diffraction modulators was successfully applied for the first time. 

The device for measuring the emission anisotropy! was also built in collaboration 
with R. Bauer (Bauer and Rozwadowski 1961). It makes possible the measurement of small 
values of the anisotropy of weak fluorescence. Values of the emission anisotropy of the 
order of 1% were determined to within 5% of the value measured. 

The relative fluorescence yields were measured by the method proposed by D. Frąc- 
kowiak and T. Marszałek (1960, 1961). 

The light source consisted throughout of a 12 V 50 W bulb supplied from a battery of 


accumulators. 


Results of the measurements, and their interpretation 


Fig. 2 shows the absorption and emission spectra of an aqueous solution of fluorescein 
in their dependence on the pH of the medium, at constant concentration of the dye. For 
the absorption spectra, the extinction coefficients are plotted on the axis of ordinates (all 
vessels had the same thickness), whereas the emission spectra are normalized to the 
same height of the maximum. No differences were found to occur in the position of the 
emission and absorption bands at the same pH values obtained by adding H,SO, 
or HNO,. However, the value of the coefficient of extinction, at identical values of pH, is 


1 Emission anisotropy (Jabłoński 1960) on excitation with polarized light is defined as follows: 


wherein 7 iI and I are the intensity components parallel and perpendicular to the electric vector of the exciting 


I 4 Ee 
light respectively, 0 — = is the degree of depolarization, and P = ee — the degree of polarization. On 
|| ISL 


the other hand, on excitation with natural light, we have 


Dar হস হন EYA 


wherein J’ and I” are the components perpendicular and parallel to the direction of the exciting light beam 
and o, and P, are the rates of depolarization and polarization at excitation with natural light, respectively. 
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slightly higher in sulfuric than in nitric acid. The temperature effect of the absorption 
spectrum of fluorescein is more marked in HNO, than in H,SO, (Fig. 5). It seems probable 
that Grigorian et al. (1960) would have obtained identical changes in the emission spectra 


A (ma) > A (mu) — 
400 450 500 550 60 400 450 500 550 600 


2 diwe n = a = rm 
450 500 550 600 0 450 500 550 600 
A(ma)-- A (mu) — 


400 


Fig. 2. Absorption and emission spectra of fluorescein solutions, for different values of pH (against the standard 
spectra at pH=10.5, shown by the dashed line): a) 12 n NH„ b) pH=7.57 (NH), c) pH—6.3 (708), d) pH 
4.15 (HNO,), e) pH=3.21(HNO,), f) pH=1.92(HNO,), g) 0.1 n H,SO,, h) 2n H,SO, i) 5n H550,, j) 10 n H,SO, 
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from sulfuric acid and nitric acid if they had applied higher H,SO, concentrations. The 
appearance of new bands in the absorption or emission spectrum points to the existence of 
different kinds of ions at different values of pH. Zanker and Peter (1958) proposed seven 
distinct possible forms of the fluorescein molecule (Fig. 3). 

It is a well-established fact that, in strongly basic solutions, bivalent negative ions are 
present (VI), whereas strongly acid solutions contain univalent positive ions (I). Hence the 
absorption band presenting 4,,,22495 mu should be attributed to modification VI, and the 


VII 
Fig. 3. Ionic modifications of fluorescein, according to Zanker. I — univalent positive ion, II and III — neutral 
molecule, IV and V — univalent negative ion, VI — bivalent negative ion, VII — modification being 


a transition between I and II 


band with Amar 27434 mu — to modification I. Modification II is colourless because of the 
oxygen bridge at the lower benzene ring. Modification VII is a transition between I and II. 
As modifications VI and V present a similar structure of the upper part, their spectra are 
undistinguishable. The spectral position of the 0,0 transitions for modifications I and VI 
(or V) is in very good agreement with the results obtained by the FEMO method?. The 
band with 4,,,4:460 mu appearing in acid solutions (from pH —7 to pH —2) should be 
attributed. to univalent negative ions (IV) or neutral molecules (III). Nevertheless, the fact 
that on adding acid to the solution an absorption band appears at first at the long-wave 


2 The calculations were carried out by Mr. A. Zieliński in the Chair of Theoretical Physics of the N. Coper- 


» 


icus University at Toruń. 
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side remains incomprehensible. It is only in very acid solutions that a band corresponding 
to cations appears to the short-wave side (the spectrum becomes similar to the one obtained 
in vitrous boric acid). Grigorian, Kantardjian and Chirkinian failed to notice that the changes. 
in the emission spectrum do not occur simultaneously with the changes in the absorption 
spectrum. This results clearly from Fig. 2. The short-wave absorption maximum 
relating to cations appears in solutions of considerably lower acidity than for emission. Even 
in 10 n H,SO, solutions emission still exhibits the maxima dominating in slightly acid and 
basic solutions, whereas already at pH —2 the long-wave absorption bands vanish completely. 

Weber (1931) observed a similar effect in 1-naphtylamin-4-sulphonate, later and Fórster 
(1950, 19508), Weller (1952, 1954, 1955, 1956, 1957), Kokubun (1958) and Derkacheva. 
(1960) for a number of other compounds. Fórster, moreover, proposed an explanation 
of the effect. He assumed that the molecule in the excited state can give off or capture 
a proton (protolytic reactions) and emit radiation as a different ion or neutral molecule. 
Such a reaction is greatly facilitated if the H* is given off by O or N. Quite generally, 
the protolytic reaction can be written as follows: 


(AH)" + B" zx A"-1+- (HBy"'1, 


wherein n and m denote the number of positive elementary charges. The energy diagram 
of the reaction is shown on Fig. 4. AE is the excitation energy of a molecule in the non- 
-dissociated form, A E' — that of the dissociated molecule, whereas Ep and Ep. are the energy 
of dissociation in the ground state and excited state, respectively. 


AE 


(AH)^ 


Fig. 4. Energy levels of molecule previous and subsequent to giving off a proton 


According to whether AE’ is smaller or larger than AE, the emission spectrum shifts. 
towards longer or shorter wavelengths. However, the effect of protolytic reactions on the 
emission spectra becomes accessible to observation only if the mean lifetime (v) of the mol- 
ecule in the excited state is longer than the mean time required by the molecule in the 
excited state for giving off the proton. As proved by investigations on proton mobility 
(Ackermann 1961, Samoylov 1957), the mean time during which a proton remains with 
a water molecule in acid media is of the order of 10-13 sec, which is by over four orders 
shorter than 7. The decisive factor here is the value of the energy of dissociation in the excited 
state which, however, should not be estimated to be much larger in fluorescein than in 
water, since even in the ground state the cations are not stable, as proved by the effect of the 
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temperature on the absorption spectra (Fig. 5). The area under the absorption curve of the 
cations diminishes strongly with increase of temperature, whereas the absorption band of the 
bivalent negative ions only undergoes a displacement towards longer wavelengths owing to the 
change in population of the vibrational levels. The value of Eh, decreases as the value of pH 
at which protolytic reactions cease to occur rises. From Fig. 2 anions are seen to exist even 
in 10 n HSO, in the excited state, which points to a considerable difference between En, 
and Ep. Although it is not possible to determine the dissociation constants for the excited 


A (mu) —— A (mu) —— 
400 420 440 460 480 500 520 400 420 440 460 480 500 520 


1:2 


400 420 440 460 480 500 520 400 420 440 460 480 500 520 
A (mu) = A (mu) —— 


Fig. 5. Effect of the temperature on the absorption spectrum: ————— coefficients of extinction at t= 15C 
— — — — coefficients of extinction at t = 40°C — . — . —. coefficients of extinction at 70°C a) pH = 11.85 


(NH,), b) pH =6.48 (HNO,), c) pH = 1.9 (H,SO,), d) pH = 1.88 (HNO,) pH was measured at £ = 20°C in the 


solutions 


state in fluorescein because of the lacton form occupying an intermediate position between 
the cations and anions, the occurrence of protolytic reactions would seem to be beyond doubt 
also in the case of fluorescein. Dissociation constants were found in simpler cases only (Weller 
1952, 1957, Derkacheva 1960). 

Figs. 6 and 7 and Tables I and II contain the values of 7 and those of the relative 
quantum yields 7, versus the pH of the medium. In the Tables I and II, the values 
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TABLE I 
NH; HNO, H50, 
pH t x 10? sec pH 1 X 10? sec pH CH,SO, rx10> sec 
11.92 4.53 6.70 | 3.60 5.33 3.40 
10.50 4.54 6.30 3.49 3.05 3.43 | 

9.90 4.52 4.62 3.43 2.13 3.21 

7.57 4.40 3.40 3.16 <1.00 3.10 

1.94 3.24 0.5 n 3.65 

<1.00 3.11 2n 3.75 

5n 3.02 

10n 4.41 

TABLE II 
NH; HNO, | H,SO, 
سسب‎ | 
| 3 7 ] | m 
pH Ãexcit Ãexcit = pH Aexcit— | Aexcit= | pH Aexcit7 Aexcit= 
| =437 mu | =490 my —437 mu | =490 ma | | —437 mu | =490 mu 

11.85 1.01 0.99 7.30 0.93 1.00 | 6.12 0.42 0.78 
11.23 1.01 0.99 6.74 0.68 096 | 602 0.40 0.74 
10.15 1.02 1.00 5.56 0.39 059 | 50 0.37 0.59 
8.50 0.98 1.00 4.95 0.32 0.35 | 3.05 0.29 0.30 
3.10 0.29 0.30 | 2.10 0.28 0.28 

1.90 0.28 — | 19 0.28 A 

<1(0.7) 0.28 — | <1(0.8) 0.28 b. 


of t and n, have been placed in separate columns for the different admixtures added to 
the initial slightly basic solution (7<pH<8) of fluorescein. We assumed n,=1 for the 
solution pH=10 and excitation in the absorption band of bivalent negative ions. 
The considerable differences in yield on excitation with wavelengths ۸-490 mu and 
A=437 mu in the case of solutions containing simultaneously anions and cations in the 
ground state is due to the fact that, on excitation in the absorption band of the cations, 
fluorescence of the cations undergoing dissociation prevails, whereas at excitation in the 
absorption band of the anions, fluorescence of the anions not participating in protolytic 
reactions is prevalent. r was measured at excitation through a BG 12 Schott filter, i.e. at 
simultaneous excitation of both bands. On the other hand, yield measurements were carried 
out on excitation through a monochromator. From either graph, the choice of the acid is 
seen to have little effect on v and n,. It is only at pH=3 that + differs somewhat for 
different acids. This is probably due to the effect of the difference in mobility of the hydrogen 
ions from HNO; and H580,. Unfortunately there was no possibility of measuring v and n, for 
each modification separately because of superposition of the emission bands, except for 
the case of bivalent negative ions. This is why the interpretation of the results becomes 


— d P 


1013 


SĘ 107 3n DEE DENS ATE De 07177777797:17111517519 
Cy, 504 
Fig. 6. Mean decay times of fluorescence versus pH of the solutions: E] — in solutions containing NH, 
© — in solutions containing HNO;, A — in solutions containing H,SO, 


OP AD RO ESCH STDC OT ORO OCHA 
pH 


Fig. 7. Relative fluorescence yield versus pH of the solutions: [5] — solutions with NH,, excitation with 

À = 490 mu, 18 — solutions with NH, excitation with A= 437 ma, © — solutions with HNO,, excitation with 

À = 490 mu, © — solutions with HNO,, excitation with A = 437 mu, A — solutions with H,SO,, excitation 
with A = 490 mu, A — solutions with H,SO,, excitation with A= 437 mu 


difficult. If, indeed, the different modifications of fluorescein possess intrinsically different 
t and 7%, the experimental results are a superposition of two independent effects, namely, 


of the intrinsically different + and 7 and of the effect of the protolytic reactions. 
In the case of 3-aminopyren-5,8,10-trisulphonate, Hanle and Scharmann (1952) found 


that when the molecule gives off a proton no change intrinsic in t and 7 occurs. Hence, on the 


3 74 denotes the absolute quantum yield, and 7, — the relative quantum yield. 
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assumption of equal r and 4 for all modifications of the fluorescein ion also, 116 6 
reactions alone suffice for explaining the observed changes in 7 and n,. Changes 
of t and 4, set in simultaneously when short-wave bands corresponding to positive and 
univalent negative ions appear in the absorption spectrum. As further changes towards lower 
pH values occur, the variations of r and n, diminish, and it is only towards higher con- 
centrations of H,SO, that t again increases, although the yield does not exhibit variations 
exceeding experimental errors (in strongly acid solutions, the error in measuring n, increases 
owing to the considerable changes in the emission spectrum). 

Let us consider what can happen to a positive or univalent negative ion after excitation. 
In a time of the order of 10-18 sec it gives off one or several protons. By giving off a single 
proton, the cation passes into a nonfluorescent lacton form. This process does not 
affect the mean decay time of the fluorescence, but solely the yield. If two or three protons 
are given off, a fluorescent negative ion arises. However, owing to reversibility of the protolytic 
reactions, there exists a certain probability (per unit time) of capture of a proton 
by the negative ion, which can lead to a non-fluorescent neutral molecule. Thus, inverse 
reactions tend to shorten the lifetime of the molecule in the excited state. However, at 
high values of the hydrogen ion concentration, a short-wave band appears in the emission 
spectrum, indicating to emission by fluorescein in the form of positive ions. The reactions, 
however, have not yet ceased completely, as the maxima corresponding to emission by 
negative ions have as yet not disappeared at the long-wave side. The lack of an appreciable 
increase in the yield also points to continued formation of the lacton form. It is noteworthy 
that the shortening of + amounts to about one third of its value for bivalent negative ions, 
whereas the yield falls to one third. This is explained by the greater probability for the for- 
mation of neutral molecules than for that of negative ions, as the latter involve giving off 
more than one proton by the cations. 

The measurements also covered the dependence of the emission anisotropy of fluo- 
rescence r, on the pH of the medium at excitation in the cation absorption band (through 
a HgM-D filter) and in that of bivalent negative ions (through an interference filter of 
Ãmax 7413 mu). The experimental results are shown in Fig. 8 and Table III. For basic solu- 


tions, r, is seen not to depend on the wavelength of the exciting light, to within experimental 


TABLE III 


NH; HNO, HSO, 
i rp(%) sj ro(96) r9(%) 
pH Aexcit = Aexcit = pH Aexcit7 Žexcit= pH Žexcit™ Aexcit ^ 
—410 mu | —473 mu —410 my | =473 mu —410 mu | =473 mu 
11.88 — 1.04 6.72 i et 1.10 6.27 _ 1:43 
10.73 1.03 1.01 6.58 1.08 1212 3.16 1.14 1.39 
10.00 — 1.02 4.36 187 1.34 2:21 1.19 — 
9.23 1.03 1.02 4.25 1.27 1.38 1.13 1:23 _ 
7.43 1.03 1.09 3.33 1.29 1.45 1.30 1111 — 
1.30 1.24 — = 
1.00 MIE — = 
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error, in accordance with the polarization spectrum for fluorescein at dye concentration 
c=10% g/cm? (Feofilov 1959). On the other hand, at pH<7, an increase in the hydrogen 
ion concentration produces an increase in the value of ry. This is due essentially to the 
shortening of +. Protolytic reactions involve a change in sign of the charge carried by the 
fluorescein ions, thus affecting the hydratation envelope. Thus, the molecule obtains a transient 
enhanced freedom of motion; this, in turn, tends to reduce ry. This effect manifests itself by 
the distinct dependence of r, on the exciting wavelength for acid solutions, and is dependent 


on the properties of the acid utilized. According to Hanle and Scharmann (1952), the process 


WY e S €. ła 4 bl e od «b W ZA KZ) 
- pH 


Fig. 8. Anisotropy of emission versus pH of solution: E] — in solution with NH; excited with A = 473 mu: 

B — in solutions with NH, excited with A = 410 mu, © — in solutions with HNO, excited with A = 473 mu 

© — in solutions with HNO, excited with A = 410 mu, A — in solutions with H,SO, excited with A = 473 my 
A — in solutions with H,SO, excited with A = 410 mu 


of giving off a proton does not affect the orientation of the molecule. The decrease in 
the degree of polarization by 3% of the measured value, as found by these authors, was 
attributed to experimental error. 

In solutions containing nitric acid, r, is maximum at pH~3, which is not the case 
with sulfuric acid. This is no doubt due to the effect of the NOS ions, which tend to enhance 
the mobility of the water molecules (Samoylov 1957, Andrussov 1958) and, at higher concen- 
trations, disrupt the hydratation envelope of the fluorescein ions. Yet in sulfuric acid 
the strongly hydratated SO” ions gradually diminish the mobility of the dye molecules. 
It is noteworthy that, at pH~3 when the greatest difference in r, occurs as between sulfuric 
and nitric acid, differences in the mean decay times are also observed to occur (Fig. 5). 
Fig. 4 shows that the effect of the temperature on the absorption spectrum is greater in nitric 
than in sulfuric acid, too. Also, the area of the absorption curve for cations is smaller in 
nitrie acid, at equal values of pH. This points to a lesser relative concentration of the fluores- 
cein cations in HNO, and, thus, toa greater number of the univalent negative ions, which are 
excited through the two filters used in the measurements of r, and contribute to the rise in r, 
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in solutions of pH<<7. This explains the fact that higher values of r„ are obtained in solutions 

with nitric acid than in those with sulfuric acid at pH ranging from 7 to 2. ¬4 
The exceedingly low value of r„ of the aqueous solution is primarily due to Brownian 

rotations and torsioned vibrations of luminescent molecules. whose effect greatly exceeds 


that of the remaining depolarizing factors. 


Results 


1. Simultaneous existence of different ions and the lacton form of fluorescein in solu- 
tion throughout a wide interval of pH values is confirmed. 

2. Protolytic reactions were found involving the loss of more than one proton by the 
excited molecule. 

3. The assumption of protolytic reactions made possible the explanation of the experi- 
mental effect of pH on the emission spectrum, mean decay time, and fluorescence yield. 

4. The effect of NO; ions on the mobility of water molecules found in direct measure- 
ments of the viscosity manifests itself also in the results obtained in measuring the emission 


anisotropy. 


5. The change in the charge of the dye ion in the excited state leads to smaller values of r,. 


The author wishes to express his great indebtedness to Professor A. Jabłoński for 
drawing his attention to the present subject and for his guidance throughout the investiga- 
tion. The author is also indebted to Dr W. Woznicki and Dr S. Poczopko for their valuable 
discussions. | 
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LETTERS TO THE EDITOR 


SOME REMARKS ON THE THEORY OF PHOTODISINTEGRATION 
OF Be AND Li? 


By A. KOWALSKA 


Institute of Theoretical Physics, Jagellonian University, Cracow 
(Received September 6, 1961) 


Some corrections in the theory of the reactions Be? (y, n) Be and Li” (y, H3) He? are 


introduced. 


In the paper of E. Guth and J. Mullin [1] on photodisintegration of Be? the fol- 
owing formula is used in order to describe the cross section for the P— D transition: 


Opp = (128/729) (n2e?/he) ho|Rppl? 


where 
R Bı )2/( (156 5 
PDT [10054 Fr) (OPH ۷۲٤٣٢ cot £r)? 
1272 FE Or) — a?r((a?4- k?) 28(1 47 ary) — COE — =] n 
X TQ (a? + k2)? (52 — 8 
^, [(a2-- k?) (1+ aro + 3(a2/k3)) + 2a2: _ (£* — p?) (L+ aro + (62/8) + 2a? 
oz E een Cd 
'242(1+ rę) OE + K?) 
| | ne En I} 


where Q = 1 — £rg cot rọ and P = 1 — 8415. 
This expression becomes infinite for &>ß. In order to obtain a correct formula changes 
must be made in two points. First P must be equal to 1 — (42/62) instead of 1 — (8412) 
second, one of the ingredients, namely 
E þ2) (1+ aro a 
(FPT 
(1019) 


= 
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must be changed into 


780 SE (1+ aro-+ (3a2/6%)) ZA 
Û (2 p 


In the new Rp the part containing (6? — 6?) in the denominator (hich forms an ex- 


0 


pression of the type E for £— Û stays finite for £-৯19 and equals: 


3(1+ arg)? (B?/a3) + 15 (1+ arp) + 5a?r6 
Ap? 


The last formula may be obtained either by putting é = f into the integral expres- 


oo 
sion Rpp = | Ri Ria, ndr, or by passing to the limit ë-> in the corrected expression 
0 


for Rpp 

These corrections change the graph on figure 3. in the cited paper [1]. The max- 
imum of the curve is at a photon energy of about 7 MeV instead of about 5 MeV 
and the maximum value of opp is ca 20 x 1075 cm?, instead of ca 15 x 107% cm? 

For the same reason the conclusions about the potential well depth in the paper 
of W. Czyż [2] on photodisintegration of Li? must be changed. The author suggests 
the existence of repulsive forces between the triton and « particle in the D state. The 
corrections given above show rather that the attracting forces are more suitable for fit- 
ting the theoretical curve to the experimental one. 
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NOTE ON POSITRON RADIATION FROM Pr” 


By S. CHOJNACKI, J. KOPYSTYŃSKI 
Institute of Experimental Physics, Warsaw University, Warsaw, 
Z. PREIBISZ, R. SOSNOWSKI, J. ZYLICZ 
Institute for Nuclear Research, Polish Academy of Sciences, Warsaw, 


I. YUTLANDOV 
Joint Institute for Nuclear Research, Dubna, USSR 


(Received June 1, 1961) 


The positron spectrum of Pr!4° was investigated with a magnetic long lens spectrometer 
It was established that the Kurie plot is a straight line from 350 keV up to the maximum energy 
of 2366-24 keV. 


In a decay scheme of Pr1% proposed by Dzhelepov et al. (1959a) there are three positron 
transitions f, a, Ps corresponding respectively to the ground state (07), to the first excited 
state (2+) with energy 1597 keV, and to the second excited state (01) with energy 1902 keV 
of 09149. As it results from measurements of the intensities of the 1597 keV and 
1902 keV transitions (Dzhelepov et al. 1959b) the components 198 and f3 are so weak 
that they cannot cause distortion in the Kurie plot of the f, transition. However, in the 
works as far published some non-linearity was observed which might suggest the existence 
of a positron component, but there is no place for such a transition in the decay scheme 
of Dzhelepov et al. (1959a). Moreover the discrepancies in the maximum energy determina- 
lions of different authors reach 10%. For these reasons it was decided to investigate the 
positron spectrum of 10140 once again. 

The measurements were carried out with a long lens spectrometer in which helical 
baffles were applied in order to separate the positron and electron radiation (Chojnacki 
et al. 1960). The investigated isotope was obtained from a neodymium fraction separated 
from a tanialum target by the chromatographic method. The target was irradiated with 
660 MeV protons (synchrocyclotron of the Joint Institute for Nuclear Research at Dubna). 
This fraction contains [৭140 which decays into Pr1%0 by electron capture. The source was 
prepared by evaporating the drops of the neodymium fraction on a thin (1.36 mg/cm?) 
aluminium foil moistened with insulin solution. The figure presents the Kurie plot for one 
of the measurements. The deviations from a straight line below 350 keV seem to be caused 


(1021) 
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by apparatus effects. The final results of this work together with the results of other اد‎ 
as well as the results predicted by semiempirical equations of Cameron (1957) and Levy 
(1957) are collected in the following table. 


Maximum energy MALONE Linearity of Kurie 
Authors of B, component | دس سه که‎ plot from Emax to 
Pons in keV | (in keV) 
1. Browne et al. (1952) 2230+20 magnetic 
spectrometer 1000 
2. Gromov et al. (1958) 2470 magnetic 
spectrometer 750 
3. Dzhelepov et al. (1959b) 2380+10 magnetic 
spectrometer 500 
4. Brabec et al. (1960) 2318+10 magnetic 
spectrometer ? 
5. Boley (1953) 2400 scintillation 
spectrometer — 
6. Gorodinsky et al. (1957) 2300 scintillation 
spectrometer — 
7. Cameron (1957) 2109 — — 
8. Levy (1957) 2810 — — 
9. Present work! 2366 +24 magnetic 
spectrometer 350 
200 
NEJ 
705 
150 
100 
50 
0 


0 300 1000 1300 2000 2500 


Kurie plot of positron radiation from 0 


1 See also Chojnacki et al. (1960). 
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